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This paper represents a detailed instruction manual for constructing the Landau expansion for 
magnetoelectric coupling in incommensurate ferroelectric magnets. The first step is to describe 
the magnetic ordering in terms of symmetry adapted coordinates which serve as complex valued 
magnetic order parameters whose transformation properties are displayed. In so doing we use the 
previously proposed technique to exploit inversion symmetry, since this symmetry had been uni- 
versally overlooked. Having order parameters of known symmetry which describe the magnetic 
ordering, we are able to construct the trilinear interaction which couples incommensurate mag- 
netic order to the uniform polarization in order to treat many of the multiferroic systems so far 
investigated. The role of this theory in comparison to microscopic models is discussed. 
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I. INTRODUCTION 

Recently there has been increasing interest in the 
interaction between magnetic and electric degrees of 
freedom, 1 Much interest has centered on a family of 
multiferroics which display a phase transition in which 
uniform ferroelectric order appears simultaneously with 
incommensurate magnetic ordering. Early examples of 
such a system whose ferroelectric behavior and mag- 
netic structure have been thoroughly studied are Ter- 
bium Manganate, TbMn0 3 (TMO)S and Nickel Vana- 
date, Ni 3 V 2 8 (NVO}^£i. A number of other systems 
have been shown to have combined magnetic and ferro- 
electric transitions^SiiSiii but the investigation of their 
magnetic structure has been less comprehensive. Initially 
this combined transition was somewhat mysterious, but 
soon a Landau expansion was developed^ to provide a 
phenomenological explanation of this phenomenon. An 
alternative picture, similar to an earlier result^ based 
on the concept of a "spin-current," and which we refer 
to as the "spiral formulation,"— has gained popularity 
due to its simplicity, but as we will discuss, the Landau 
theory is more universally applicable and has a number 
of advantages. The purpose of the present paper is to 
describe the Landau formulation in the simplest possi- 
ble terms and to apply it to a large number of currently 
studied multiferroics. In this way we hope to demystify 
this formulation. 

It should be noted that this phenomenon (which we 
call "magnetically induced ferroelectricity" ) is closely 
related to the similar behavior of so-called "improper 
ferroelectrics," which are commonly understood to be 
the analogous systems in which uniform magnetic 
order (ferromagnetism or antiferromagnetism) drive 
ferroelectricity— Several decades ago such systems were 
studied^ and reviewecU*i£ and present many parallels 
with the recent developments. 

One of the problems one encounters at the outset is 
how to properly describe the magnetic structure of sys- 
tems with complicated unit cells. This, of course, is a 
very old subject, but surprisingly, as will be documented 



below, the full ramifications of symmetry are not widely 
known. Accordingly, we feel it necessary to repeat the 
description of the symmetry analysis of magnetic struc- 
tures. While the first part of this symmetry analysis 
is well known to experts, we review it here, especially 
because our approach is often far simpler and less tech- 
nical than the standard one. However, either approach 
lays the groundwork for incorporating the effects of in- 
version symmetry, which seem to have been overlooked 
until our analysis of NVQ 3 i 4 i 5 i 6 i 7 and TMO. 3 Inversion 
symmetry was also addressed by Schweizer with a subse- 
quent correction— Very recently a more formal approach 
to this problem has been given by Radaelli and Chapon. 18 
But, at least in the simplest cases, the approach initially 
proposed by us and used here seems easiest. We apply 
this formalism to a number of currently studied multi- 
ferroics, such as MnW0 4 (MWO), TbMn 2 5 (TM025), 
YMn 2 5 (YM025), and CuFe0 2 (CFO). As was the case 
for NVQ 4 i 5 i 6 i 7 and TMO) 3 - once one has in hand the sym- 
metry properties of the magnetic order parameters, one 
is then able to construct the trilinear magnetoelectric 
coupling term in the free energy which provides a phe- 
nomenological explanation of the combined magnetic and 
ferroelectric phase transition. 

This paper is organized in conformity with the above 
plan. In Sec. II we review a simplified version of the sym- 
metry analysis known as representation theory, in which 
we directly analyze the symmetry of the inverse suscep- 
tibility matrix. Here we also review the technique we 
proposed some time ago&i^iSi 7 . to incorporate the con- 
sequences of inversion symmetry. In Sec. Ill we apply 
this formalism to develop magnetic order parameters for 
a number of multiferroic systems and in Eq. i|122|l we 
give a simple example to show how inversion symme- 
try influences the symmetry of the allowed spin distribu- 
tion. Then in Sec. IV, we use the symmetry of the order 
parameters to construct a magnetoelectric coupling free 
energy, whose symmetry properties are manifested. In 
Sec. V we summarize the results of these calculations 
and discuss their relation to calculations based on the 
spin current modeU^ or the phenomenology of contin- 
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uum theory^ 



II. REVIEW OF REPRESENTATION THEORY 

As we shall see, to understand the phenomenology of 
the magnetoelectric coupling which gives rise to the com- 
bined magnetic and ferroelectric phase transition, it is es- 
sential to characterize and properly understand the sym- 
metry of the magnetic ordering. In addition, as we shall 
see, to fully include symmetry restrictions on possible 
magnetic structures that can be accessed via a continu- 
ous phase transition is an extremely powerful aid in the 
magnetic structure analysis, Accordingly in this section 
we review how symmetry considerations restrict the pos- 
sible magnetic structures which can appear at an order- 
ing transition. The full symmetry analysis has previ- 
ously been presented elsewhere but it is useful 
to repeat it here both to fix the notation and to give 
the reader convenient access to this analysis which is so 
essential to the present discussion. To avoid the com- 
plexities of the most general form of this analysis (called 
representation theory) ^Li&iS we will limit discussion to 
systems having some crucial simplifying features. First, 
we limit consideration to systems in which the magnetic 
ordering is incommensurate. In the examples we choose k 
will usually lie along a symmetry direction of the crystal. 
Second, we only consider systems which have a center of 
inversion symmetry, because it is only such systems that 
have a sharp phase transition at which long-range ferro- 
electric order appears. Thirdly, we restrict attention to 
crystals having relatively simple symmetry. (What this 
means is that except for our discussion of TbM^Os we 
will consider systems where we do not need the full ap- 
paratus of group theory, but can get away with simply 
labeling the spin functions which describe magnetic order 
by their eigenvalue under various symmetry operations.) 
By avoiding the complexities of the most general situ- 
ations, it is hoped that this paper will be accessible to 
more readers. Finally, as we will see, it is crucial that 
the phase transitions we analyze are either continuous 
or very nearly so. In many of the examples we discuss, 
our simple approach^ is vastly simpler than that of stan- 
dard representation theory augmented by representation 
theory specialized techniques to explicitly include exploit 
inversion symmetry 



A. Symmetry Analysis of the Magnetic Free 
Energy 

In this subsection we give a review of the formalism 
used previously^ and presented in detail in Refs. 
Since we are mainly interested in symmetry properties, 
we will describe the magnetic ordering by a version of 
mean-field theory in which one writes the magnetic free 



energy F M as 

r,a;r'/3 

+C(5 4 ) , (1) 

where S a (r) is the thermally averaged a-component of 
the spin at position r. In a moment, we will give an 
explicit approximation for the inverse susceptibility x- 
We now introduce Fourier transforms in either of two 
equivalent formulations. In the first (which we refer to 
as "actual position" ) one writes 

S a (q,r) = AT- 1 ^5 a (R + r)e^( R + T ) (2) 

R 

whereas in the second (which we refer to as "unit cell") 
one writes 

S Q (q,r) = N-^S^R + t)^^ , (3) 

R 

where N is the number of unit cells in the system, r 
is the location of the rth site within the unit cell, and 
R is a lattice vector. Note that in Eq. (0 the phase 
factor in the Fourier transform is defined in terms of the 
actual position of the spin rather than in terms of the 
origin of the unit cell, as is done in Eq. J2Jl. In some 
case (viz. NVO) the results are simpler in the actual 
position formulation whereas for others (viz. MWO) the 
unit cell formulation is simpler. We will use whichever 
formulation is simpler. In either case the fact that S a 
has to be real indicates that 

S a (-q,r) = S a (q,r)* ■ (4) 

We thus have 

Fm = \ X^(q;T,r')^(q,r)*^(q,r') 

q;r,r',Q./3 

+C(5 4 ) , (5) 
where (for the "actual position" formulation) 

X^(q;r,r') = £ X ^(r, R + /^-(R+t'-t) (g) 

R 

To make our discussion more concrete we cite the sim- 
plest approximation for a system with general anisotropic 
exchange coupling, so that the Hamiltonian is 

H = ^/3(r,r')sa(r) S/3 (r')+^A^ SQ (r) 2 ,(7) 

Q,/3;r,r' ar 

where s a (r) is the a-component of the spin operator at 
r and we have included a single ion anisotropy energy 
assuming three inequivalent axes, so that the K a are all 
different. One has that 

*J(r,r') = J a p(r,r') + [K a + ckT]5 a ,p5 r , T , , (8) 
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where 8 a ^ is unity if a — b and is zero otherwise and c 
is a spin-dependent constant of order unity, so that ckT 
is the entropy associated with a spin S. Then if we have 
one spin per unit cell, one has 

Xa/Kq) = 5 a p^2Ji[cos(a a q x ) + cos(a a q y ) 

+ cos(a a q z )} + akT + K a ^j , (9) 

where a a is the lattice constant in the a-direction2S and 
we assume that K x < K y < K z . Graphs of x _1 (q) are 
shown in Fig. 1 for both the ferromagnetic ( J\ < 0) and 
antifcrromagnetic ( J\ > 0) cases. For the ferromagnetic 
case we now introduce a competing antiferromagnetic 
next-nearest neighbor (nnn) interaction J2 > along the 
x-axis, so that 

Xaailx, % = 0, q z = 0) = [4Ji + 2Ji cos(a x q x ) 
+2J 2 cos{2a x q x ) + akT + K a ] , (10) 

and this is also shown in Fig. 1. As T is lowered one 
reaches a critical temperature where one of the eigen- 
values of the inverse susceptibility matrix becomes zero. 
This indicates that the paramagnetic phase is unstable 
with respect to order corresponding to the critical eigen- 
vector associated with the zero eigenvalue. For the fer- 
romagnet this happens for zero wavevector and for the 
antiferromagnet for a zone boundary wavevector in agree- 
ment with our obvious expectation. For competing in- 
teractions we see that the values of the J's determine 
a wavevector at which an eigenvalue of is minimal. 
This is the phenomenon called "wavevector selection," 
and in this case the selected value of q is determined by 
extremizing % _1 to be2i 

cos(a x q) = -Ji/(4J 2 ) , (11) 

providing J2 > — Ji/4. (Otherwise the system is ferro- 
magnetic.) Note also, that crystal symmetry may select 
a set of symmetry-related wavevectors, which comprise 
what is known as the star of q. (For instance, if the sys- 
tem were tetragonal, then crystal symmetry would imply 
that one has the same nnn interactions along the y-axis, 
in which case the system selects a wavevector along the 
x-axis and one of equal magnitude along the y-axis. 

From the above discussion it should be clear that if we 
assume a continuous transition so that the transition is 
associated with the instability in the terms in the free 
energy quadratic in the spin amplitudes, then the nature 
of the ordered phase is determined by the critical eigen- 
vector of the inverse susceptibility, i. e. the eigenvector 
associated with the eigenvalue of inverse susceptibility 
which first goes to zero as the temperature is reduced. 
Accordingly, the aim of this paper analyze how crystal 
symmetry affects the possible forms of the critical eigen- 
vector. 

When the unit cell contains n > 1 spins, the inverse 
susceptibility for each wavevector q is a 3n x 3n matrix. 



The ordering transition occurs when, for some selected 
wavevector (s), an eigenvalue first becomes zero as the 
temperature is reduced. In the above simple examples 
involving isotropic exchange interactions, the inverse sus- 
ceptibility was 3x3 diagonal matrix, so that each eigen- 
vector trivially has only one nonzero component. The 
critical eigenvector has spin oriented along the easiest 
axis, i. e. the one for which K a is minimal. In the 
present more general case n > 1 and arbitrary interac- 
tions consistent with crystal symmetry are allowed. To 
avoid the technicalities of group theory, we use as our 
guiding principle the fact that the free energy, being an 
expansion in powers of the magnetizations relative to the 
the paramagnetic state, must be invariant under all the 
symmetry operations of the crystallise This is the same 
principle that one uses in discussing the symmetry of the 
electrostatic potential in a crystal^ We now focus our 
attention on the critically selected wavevector q which 
has an eigenvalue which first becomes zero as the tem- 
perature is lowered. This value of q is determined by 
the interactions and we will consider it to be an experi- 
mentally determined parameter. Operations which leave 
the quadratic free energy invariant must leave invariant 
the term in the free energy F2 (q) which involves only the 
selected wavevector q, namely 

^(q) = \ E ^(q;r,r')5 a (q ) r)*5 /3 (q ) T')(.12) 

r,r' ,a,/3 

Any symmetry operation takes the original variables be- 
fore transformation, S a (q, r) , into new ones indicated by 
primes. We write this transformation as 

^(q,r) =J2 U <*T; a 'T>S a >(<i,T') . (13) 

a'T 

According to a well known statement of elementary 
quantum mechanics, if an operator T commutes with 
X _1 (q), then the eigenvectors of X _1 (q) are simulta- 
neously eigenvectors of T. (This much involves a well 
known analysis ^SLSl We will later consider the effect 
of inversion, the analysis of which is universally over- 
looked). We will apply this simple condition to a num- 
ber of multiferroic systems currently under investigation. 
(This approach can be much more straightforward than 
the standard one when the operations which conserve 
wavevector unavoidably involve translations.) As a first 
example we consider the case of NVO and use the "ac- 
tual position" Fourier transforms. In Table U we give 
the general positions (this set of positions is the so-called 
Wyckoff orbit) for the space group Cmca (#64 in Ref. 
l28l) of NVO and this table defines the operations of the 
space group of Cmca. In Table HTl we list the positions of 
the two types of sites occupied by the magnetic (Ni) ions, 
which are called "spine" and "cross-tie" sites in recogni- 
tion of their distinctive coordination in the lattice, as 
can be seen from Fig. [21 where we show the conven- 
tional unit cell of NVO. Experiments^*^ indicate that as 
the temperature is lowered, the system first develops in- 



FIG. 1: Inverse susceptibility x~ (° 5 0, 0). a) Ferromagnetic model (Ji < 0), b) Antiferromagnetic model (Ji < 0), and c) 
Model with competing interactions (the nn interaction is antiferromagnetic). In each panel one sees three groups of curves. 
Each group consists of the three curves for x^a{q) which depend on the component label a due to the anisotropy. The x axis 
is the easiest axis and the z axis is the hardest. (If the system is orthorhombic the three axes must all be inequivalent. The 
solid curves are for the highest temperature, the dashed curves are for an intermediate temperature, and the dash-dot curves 
are for T = T c , the critical temperature for magnetic ordering. Panel c) illustrates the nontrivial wavevector selection which 
occurs when one has competing interactions. 



Er = (x,y,z) 2 c r = (x, y + 1/2, z + 1/2) 

2 b r = (x,y + 1/2, z + 1/2) 2 a r = (x,y,z) 

lv=(x,y,z) m c r = (a;, y + l/2,z + 1/2) 
m b r = (x,y + 1/2, z + 1/2) m a r = (x, y, z) 



TABLE I: General positions— *— within the primitive unit cell 
for Cmca which describe the symmetry operations— of this 
space group. 2 a is a two-fold rotation (or screw) axis and m a 
is a mirror (or glide) which takes r a into — r a . 





= (0.25,-0.13,0.25) 


r S 2 


= (0.25,0.13,0.75) 


r S 3 


= (0.75,0.13,0.75) 


r a 4 


= (0.75,-0.13,0.25) 


r c i 


(0,0,0) 


r C 2 


(0.5,0,0.5) 



TABLE II: Positions**** of Ni 2+ carrying S=l within the 
primitive unit cell illustrated in Fig. Here r sn denotes the 
position of the nth spine site and r cn that of the nth cross-tie 
site. NVO orders in space group Cmca, so there are six more 
atoms in the conventional orthorhombic unit cell which are 
obtained by a translation through (0.5a, 0.5&, 0). 



commensurate order with q along the a-direction with 
q « 0.28^4 

The group of operations which conserve wavevector are 
generated by a) the two-fold rotation 2 X and b) the glide 
operation m z , both of which are defined in Table [Q We 
now discuss how the Fourier spin components transform 
under various symmetry operations. Here primed quan- 
tities denote the value of the quantity after transforma- 
tion. Let O = O s O r be a symmetry operation which 
we decompose into operations on the spin O s and on the 
position O r . The effect of transforming a spin by such 



an operator is to replace the spin at the "final" position 
R/ by the transformed spin which initially was at the 
position C~ 1 R/. So we write 

S' a (Rf,T f ) - CSaiO-^Tf]) 

= UOs)S a CR l ,T l ) , (14) 

where the subscripts "i" and "f ' denote initial and final 
values and £, a (O s ) is the factor introduced by O s for a 
pseudovector, namely 

U2 X ) = 1. €»(2«)=C»(2 1B )=-1, 
Cx(m z ) = £ y (m z ) = -1 , £z(m 2 ) = 1 • (15) 

Note that OS a (R, r) is not the result of applying O to 
move and reorient the spin at R + r, but instead is the 
value of the spin at R + r after the spin distribution 
is acted upon by O. Thus, for actual position Fourier 
transforms we have 

S'M^f) - iV- 1 ^^(R / ,r / )e^( R /+' r /) 

R 

R 

= UOs)S a (q,T l )e^ R f+ T f-^- T ^ .(16) 

We may write this as 

OSafaTf) = ea(a)^(q,T 4 )e Jq -^ +T ^ R '- x '](17) 

This formulation may not be totally intuitive, because 
one is tempted to regard the operation O acting on a 
spin at an initial location and taking it (and perhaps re- 
orienting it) to another location. Here, instead, we con- 
sider the spin distribution. The transformed distribution 
at a location r is related to the distribution at the initial 
location O r T 1 r. 

Similarly, the result for unit cell Fourier transforms is 

S'M,T f ) = £ a (a)^(q,T;) e ^ R /- R «] . (18) 
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As before, we may write this as 

OS a (q,r f ) = UOs)S a (<i,T i y^ R i~ R 'K (19) 
Under transformation by inversion, = 1 and 

S'M,T f y = 7V- 1 ^5 ct (R,,x 4 )e-^( R /+ T ^ 



R 



= S a (q,T i )e^- R /- T '- R *- T '] 

= S a (q,n) (20) 

for actual position Fourier transforms. For unit cell 
transforms we get 

S'M,T f y = S a (q,T i )e*»-[- R '- R '] 

= S a (q,T % )e^ T f +T ^ . (21) 



Irrep 


Ti 


r 2 


r 3 


r 4 


A(2„) = 


+ 1 


+i 


-l 


-l 


A(m z ) = 


+ 1 


-l 


-l 


+i 


S(q,sl) 


n> 


n< s 
"a 
n s 


n\ 


n» 


S(q, s2) 




n< s 

~< 


< 


— 
"a 






~ n s 






S(q, S3) 


— 

"5 


n< s 

-< 


n\ 


-< 




— ri c s 


n s 


— ri c s 


< 


S(q,s4) 


—n't 

-n h s 


n< s 
"a 


n s 
n» 3 


-n» 




n s 


~ n s 






S(q,cl) 


n c 





n c 







< 





S(q,c2) 






n c 







< 
-< 




-< 



TABLE III: Allowed spin functions (i. e. actual position 
Fourier coefficients) within the unit cell of NVO for wavevec- 
tor (q, 0, 0) which are eigenvectors of 2 X and m z with the 
eigenvalues A listed. Irrep stands for the irreducible represen- 
tation as labeled in Ref. 0. The labeling of the sites is as in 
Table El and Fig. H Here n£ (p =s or c, a = a, b, c) denotes 
the complex quantity rip(q). 

Now we apply this formalism to find the actual position 
Fourier coefficients which are eigenf unctions of the two 
operators 2 X and m z . In so doing note the simplicity of 
Eq. I|16(l : since, for NVO, the operations 2 X and m z do 
not change the a; coordinate, we simply have 

S'M,r f ) = us'M-ri) ■ (22) 

Thus the eigenvalue conditions for 2 X acting on the spine 
sites (#l-#4) are 

S a (q,l)' = £ a (2 x )S a (q,2) = \{2 x )S a (q,l) 

S a {q,2)' = £ a (2 x )S a (q,l) =\{2 x )S a {q, 2) 

S a (q,3)' = ea(2x)S a (g,4)=A(2 x )5 a ( 9) 3) 

S a (q,4)' = t; a (2 x )S a (q,3) = X(2 x )S a (q, 4) , (23) 



from which we see that X(2 X ) — ±1 and 

S a (q,2) = K„(2 x )/A(2 B )]5„(g,l) 
S a (q,3) = K„(2 x )/A(2 B )]5„(g,4) 



(24) 



The eigenvalue conditions for m z acting on the spine sites 
are 

S a (q,l)' = £ a (m z )S a (q,4) = X(m z )S a (q, 1) 
S a (q,A)' = £ a (m z )S a (q, 1) = X(m z )S a (q,A) 
S a (q,2)' = £ a (m z )S a (q, 3) = X(m z )S a (q, 2) 
S a (q,3)' = t a (m z )S a (q,2) = X(m z )S a (q,3) , (25) 



from which we see that X(m z ) = ±1 and 

S a (q,4) = [€a(m z )/X(m z )]S a (q, 1) . 



(26) 



We thereby construct the eigenvectors for the spine sites 
as given m Table HTU The results for the cross-tie sites 
are obtained in the same way and are also given in the 
table. Each set of eigenvalues corresponds, in techni- 
cal terms, to a single irreducible representation (irrep). 
Since each operator can have either of two eigenvalues, 
we have four irreps to consider. Note that these spin 
functions, since they are actually Fourier coefficients, are 
complex- valued quantities. [The spin itself is real be- 
cause F(— q) = F(q)*.] Note that each column of Table 
nil gives the most general form of an allowed eigenvec- 
tor for which one has 4 (or 5, depending on the irrep) 
independent complex constants. To further illustrate the 
meaning of this table we explicitly write, in Eq. 145J1 . 
below, the spin distribution arising from one irrep, Tzt. 

So far, the present enalysis reproduces the standard 
results and indeed computer programs exist to construct 
such tables. But for multiferroics it may be quicker to 
obtain and understand how to construct the possible spin 
functions by hand rather than to understand how to use 
the program! Usually these programs give the results 
in terms of unit cell Fourier transforms, which we claim 
are not as natural a representation in cases like NVO. In 
terms of unit cell Fourier transforms the eigenvalue con- 
ditions for 2 X acting on the spine sites (#l-#4) are the 
same as Eq. l|2l^l for actual position Fourier transforms 
because the operation 2 X does not change the unit cell. 
However, for the glide operation m z this is not the case. 
If we start from site #1 or site #2 the translation along 
the y axis takes the spin to a final unit cell displaced 
by (—a/2)i + (b/2)j, whereas if we start from site #3 or 
site #4 the translation along the y axis takes the spin 
to a final unit cell displaced by (a/2)i + (b/2)j. Now 
the eigenvalue conditions for m z acting on the spine sites 
(#l-#4) are 



S a (q,l)' 
S a (qAY 
S a (q,2)' 
S a (q,3)' 



£ a (m z )S a (q,4)r) 
£a{m z )S a (q, l)r)* 
£ a {m z )S a (q,3)r) 
£ a (m z )S a (q,2)r)* 



X(m z )S a (q, 1) 

= X(m z )S a (q,4) 
X(m z )S a (q,2) 
= X(m z )S a (q,3) ,(27) 
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FIG. 2: (Color online). Ni sites in the conventional unit cell of 
NVO. The primitive translation vectors v n are vi = (a/2)a + 
(6/2)6, v 2 = (a/2)o - (6/2)6, and v 3 = cc. The "cross-tie" 
sites (on-line=blue) cl and c2 lie in a plane with 6 = 0. The 
"spine" sites (on-line=red) are labeled si, s2, s3, and s4 and 
they may be visualized as forming chains parallel to the a- 
axis. These chains are in the buckled plane with 6 = ±5, 
where 8 = 0.136 as is indicated. Cross-tie sites in adjacent 
planes (displaced by (±6/2)6) are indicated by open circles. 
Spine sites in adjacent planes are located directly above (or 
below) the sites in the plane shown. In the incommensurate 
phases the wavevector describing magnetic ordering lies along 
the a axis. The axis of the two-fold rotation about the x-axis 
is shown. The glide plane is indicated by the mirror plane at 
z = | and the arrow above m z indicates that a translation of 
6/2 in the y-direction is involved. 



where r; = cxp(z7rg). One finds that all entries for 
5(q, s3), 5(q, s4), and 5(q, c2) now carry the phase fac- 
tor r\* — exp(— iirq). But this is just the factor to make 
the unit cell result 

S(R,r) = S(q,r)e-^ R (28) 

be the same (to within an overall phase factor) as the 
actual position result 

S(R,t) = S(q,r)e-' !q -( R +' r ) . (29) 

We should emphasize that in such a simple case as NVO, 
it is actually not necessary to invoke any group theo- 
retical concepts to arrive at the results of Table IHII for 
the most general spin distribution consistent with crystal 
symmetry. 

More importantly, it is not commonly 
understood^2Si2i that one can also extract infor- 
mation using the symmetry of an operation (inversion) 



which need does not conserve wavevector 
Since what we are about to say is unfamiliar, we start 
from first principles. The quadratic free energy may be 
written as 

F 2 = E E KpSa^rrSp^y) , (30) 

q r,r';a/3 

where we restrict the sum over wavevectors to the star of 
the wavevector of interest. One term of this sum is 

F 2 (q ) = E F a^S a (qo,T)*S (q_o,T') . (31) 

r,r';a/3 

It should be clear that the quadratic free energy, F 2 is in- 
variant under all the symmetry operations of the param- 
agnetic space group (i. e. what one calls the space group 
of the crystal) i22i2^ For centrosymmetric crystals there 
are three classes of such symmetry operations. The first 
class consists of those operations which leave go invariant 
and these are the symmetries taken into account in the 
usual formulationi2£i2§i2! The second class consists of op- 
erations which take go into another wavevector of the star 
(call it qi), where qi 7^ — qo- Use of these symmetries al- 
lows one to completely characterize the wavefunction at 
wavevector qi in terms of the wavefunction for qo . These 
relations are needed if one is to discuss the possibility of 
simultaneously condensing more than one wavevector in 
the star of qiS Finally, the third class consists of spatial 
inversion (unless the wavevector and its negative differ 
by a reciprocal lattice vector, in which case inversion be- 
longs in class #1). The role of inversion symmetry is 
almost universally overlooked ^S&SL as is evident from 
examination of a number of recent papers. Unlike the 
operations of class #1 which takes S n (q) into an SVi'(q) 
(for irreps of dimension one which is true for most cases 
considered in this paper), inversion takes 5„(q) into an 
S n '(—q). Nevertheless it does take the form written in 
Eq. (|31|l into itself and restricts the possible form of the 
wavefunctions. So we now consider the consequences of 
invariance of F 2 under inversion^^iSii For this purpose 
we write Eq. I|12|) in terms of the spin coordinates m of 
Table ITTT1 (The result will, of course, depend on which 
irrep T we consider.) In any case, the part of F 2 which 
depends on qo can be written as 

F 2 (q) = E ^'^(qo,r)*^(q ,r') 

r,r';a/3 

= E G N , a ., N ,, p [ln a N ]*[ln%,] , (32) 

where N and N' assume the values "s" for spin and "c" 
for cross-tie and a and (3 label components. Now we need 
to understand the effect of I on the spin Fourier coeffi- 
cients listed in Table IIIII Since we use actual position 
Fourier coefficients, we apply Eq. (|20|) . For the cross-tie 
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variables (which sit at a center of inversion symmetry) in- 
version takes the spin coordinates of one spine sublattice 
into the complex conjugate of itself: 

ZS(q,cn) = [S(q,cn)]* . (33) 

Thus in terms of the n's this gives 

Xn a c = [n°]* , a = x,y,z . (34) 

The effect of inversion on the spine variables again fol- 
lows from Eq. I|20[) . Since inversion interchanges sublat- 
tice #1 and #3, we have 

[S(q,s3)]'= [S(q,sl)]* . (35) 

For X(2 X ) = \(m z ) = +1 (i. e. for irrep Ti), we substi- 
tute the values of the spin vectors from the first column 
of Table EH to get 

A~<] = [<]* , A<\ = [<T , 

= [<]* . (36) 

Note that some components introduce a factor — 1 under 
inversion and others do not. (Which ones have the minus 
signs depends on which irrep we consider.) If we make a 
change of variable by replacing n" in column #1 of Table 
IIIII by in™ for those components for which I introduces 
a minus sign and leave the other components alone, then 
we may rewrite the first column of Table IIIII in the form 
given in Table llVl In terms of these new variables one 
has 

I[nf\ = [<]* . (37) 

(It is convenient to define the spin Fourier coefficients 
so that they all transform in the same way under inver- 
sion. Otherwise one would have to keep track of variables 
which transform with a plus sign and those which trans- 
form with a minus sign.) Repeating this process for all 
the other irreps we write the possible spin functions as 
those of Table IIVI We give an explicit formula for the 
spin distribution for one irrep in Eq. (|45|) below. 

Now we implement Eq. i|32|) , where the spin functions 
are taken to be the variables listed in Table IIVI First 
note that the matrix G in Eq. (|32|l has to be Hermitian 
to ensure that F2 be real: 

Gm,<x;N,P = [GN,P;M,a]* ■ (38) 

Then, using Eq. we find that Eq. (J32J) is 

i^(qo) = X] \ n M]* G M,a-N,l3n P N 

= ^2 \ In M]* G M,a-NA ln N\ 

M,a;N,0 

= ^ n MGM,a;N,fi[Tlff]* 
M,a;N,0 

= E [<iYGn, P -mA^n\ > (39) 

M,a;N,0 



Irrep— 


Ti 


r 2 


r 3 


r 4 


A(2 X ) = 


+ 1 


+i 


-l 


-l 


X(m z ) = 


+ 1 


-l 


-l 


+i 




in a s 








S(q,sl) 


n\ 


■ b 

m s 


n". 


■ b 

vn s 




in c s 


n c „ 


in c s 


n s 




in" 




—iris 


—Tig 


S(q,s2) 


-W. 


-in* 


n" s 


in b s 




— in c s 




in c s 


n s 




— 


n s 


—in^ 


n s 


S(q, s3) 


< 


■ b 

— m s 


"a 






— in^ 


n c B 


—in% 


n s 




—in a s 






-Tig 


S(q,s4) 


-< 


in b s 


< 


-in\ 






~< 


—in c s 






< 


n c 








S(q,cl) 








"c 


< 










n c 








n c 








S(q,c2) 








< 


-< 










-< 


< 



TABLE IV: As Table E3] except that now the effect of in- 
version symmetry is taken into account, as a result of which, 
apart from an overall phase factor all the n's in this table can 
be taken to be real-valued. 

where, in the last line, we interchanged the roles of the 
dummy indices M, a and N, (3. By comparing the first 
and last lines, one sees that the matrix G is symmetric. 
Since this matrix is also Hermitian, all its elements must 
be real valued. Thus all its eigenvectors can be taken 
to have only real-valued components. But the m's are 
allowed to be complex valued. So, the conclusion is that 
for each irrep, we may write 

n%(T) = e^[r%(r)} , (40) 

where the r's are all real valued and 4t is an overall 
phase which can be chosen arbitrarily for each T. It is 
likely that the phase will be fixed by high-order Umklapp 
terms in the free energy, but the effects of such phase 
locking may be beyond the range of experiments. 

It is worth noting how these results should be (and in 
a few case a" ' 1 ^ have been) used in the structure determi- 
nations. One should choose the best fit to the diffraction 
data using, in turn, each irrep (i. e. each set of eigen- 
values of 2 X and m z ). Within each representation one 
parametrizes the spin structure by choosing the Fourier 
coefficients as in the relevant column of Table IIVI Note 
that instead of having 4 or 5 complex coefficients to de- 
scribe the six sites within the unit cell (see Table Mill , 
one has only 4 or 5 (depending on the representation) 
real-valued coefficients to determine. The relative phases 
of the complex coefficients have all been fixed by invok- 
ing inversion symmetry. This is clearly a significant step 
in increasing the precision of the determination of the 
magnetic structure from experimental data. 
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B. Order Parameters 

We now review how the above symmetry classification 
influences the introduction of order parameters which al- 
low the construction of Landau expansions!^ The form 
of the order parameter should be such that it has the po- 
tential to describe all ordering which are allowed by the 
quadratic free energy F^- Thus, for an isotropic Heisen- 
berg model on a cubic lattice, the order parameter has 
three components (i. e. it involves a three dimensional 
irrep) because although the fourth order terms will re- 
strict order to occur only along certain directions, as far 
as the quadratic terms are concerned, all directions are 
equivalent. The analogy here is that the overall phase 
of the spin function 6(T) is not fixed by the quadratic 
free energy and accordingly the order parameter must 
be a complex variable which includes such a phase. One 
also recognizes that although the amplitude of the critical 
eigenvector is not fixed by the quadratic terms in the free 
energy, the ratios of its components are fixed by the spe- 
cific form of the inverse susceptibility matrix. Although 
we do not wish to discuss the explicit form of this matrix, 
what should be clear is that the components of the spins 
which order must be proportional to the components of 
the critical eigenvector. The actual amplitude of the spin 
ordering is determined by the competition between the 
quadratic and fourth order terms in the free energy. If T p 
is the irrep which is critical, then just below the ordering 
temperature we write 



nft(q) = <r p (qn(r p ) 



v) i 



(41) 



where the r's are real components of the critical eigenvec- 
tor (coming from irrep T p ) of the matrix G of Eq. 1M2I) 
and are now normalized by 



aN 



(42) 



Here the order parameter for irrep T(q), <r p (q) is a com- 
plex variable, since it has to incorporate the arbitrary 
complex phase <p p associated with irrep T p : 



o>(±M) = °> e 



(43) 



The order parameter transforms as indicated in the tables 
by its listed eigenvalues under the symmetry operations 
2 X and m z . Since the components of the critical eigenvec- 
tor are dominantly determined by the quadratic terms^ 
one can say that just below the ordering temperature the 
description in terms of an order parameter continues to 
hold but 



IT,. 



(44) 

1/2 but corrections 



where mean-field theory gives j3 - 
due to fluctuation are expected^ 

To summarize and illustrate the use of Table IIVI we 
write an explicit expression for the magnetizations of the 



#1 spine sublattice and the #1 cross-tie sublattice for 
irrep T 4 [A^) = — 1 and A(m 2 ) = +1]. We use the 
definition of the order parameter and sum over both signs 
of the wavevector to get 





v si ^ 

v ) / 


— 9(T/l7 ,a ' fOSlHT -r- 


6 a) 


q, 


v si J 

v ) / 


— 9iT i SI Tl ( (IT A- i 


1>a) 


q 

u z 


v si 1 

v ) / 


— ^STaV 2 CCt^lflT -\- 

— — \i zj. i „ ^uo^y j/ i 




q 

u x 


'r s2"l 


— — 9fT/i r x COSI (IT - 


r v^4 


s, 

Uy 


r s2) 


= 2cT4/*^ SIH^QX -\- ( 


i>A ) 

J 4 J 


J Z 




— Oct a t z pn^f /i t -4- 


! J 4 ) 




v S3 1 ) 


— ^(T A T X POsfoT -I- 


6 a) 


S. 

Uy 


Y s3) 


= — 2fT/i sin ( (it - 


- 6a 


q 


I , hO) 


— rr i t ct\q.( nnr 1 


H) 






— — 9*T/i r x rrmi cit - 


\- 6a 


q 

Jy 


f r s4^1 


— 0/T i T^ Q1T1I /7T - 




Sz 


[r, s4) 


= 2gat z s cos{qx + 


h) 


S x 


(r,cl) 


= 




Sy 


(r,cl) 


= 2g^t v c cos(qx + 


fa) 


S z 


(r, cl) 


= 2u^r z c cos(qa; + 


fa) 


S x 


(r, cl) 


= 




Sy 


(r,c2) 


= — 2ct4^ cos(ga; - 


>rfa 


S z 


(r,c2) 


= 2a 4 r z cos(qx + 


fa) 



(45) 

and similarly for the other irreps. Here r = (x,y,z) is 
the actual location of the spin. Using explicit expressions 
like the above (or more directly from Table Hv|) . one can 
verify that the order parameters have the transformation 
properties: 



2 x (Ti(q) = 

2icr 2 (q) = 

2icr 3 (q) = 

2 x er 4 (q) = 



-o-i(q) , 

-CT 2 (q) , 

-CT 3 (q) , 

-CT 4 (q) , 



^cri(q) 
m z a 2 {c\) 

m z <r 4 (q) 



and 



Jcr n (q) = [<r n (q)Y 



+o-i(q) , 

-o-2(q) , 
-o- 3 (q) , 
+<T 4 (q) (46) 

(47) 



Note that even when more than a single irrep is present, 
the introduction of order parameters, as done here, pro- 
vides a framework within which one can represent the 
spin distribution as a linear combination of distributions 
each having a characteristic symmetry, as expressed by 
Eq. (|4^|l . When the structure of the unit cell is ignored^ 
that information is not readily accessible. Also note that 
the phase of each irrep T n is defined so that when <j) n = 0, 
the wave is inversion-symmetric about r = 0. For a single 
irrep this specification is not important. However, when 
one has two irreps, then inversion symmetry is only main- 
tained if their phases are equal. 

In many systems, the initial incommensurate order 
that first occurs as the temperature is lowered becomes 
unstable as the temperature is further lowered^ Typi- 
cally, the initial order involves spins oriented along their 
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easy axis with sinusoidally varying magnitude. However, 
the fourth order terms in the Landau expansion favor 
fixed length spins. As the temperature is lowered the 
fixed length constraint becomes progressively more im- 
portant and at a second, lower, critical temperature a 
transition occurs in which transverse components become 
nonzero. Although the situation is more complicated 
when there are several spins per unit cell, the result is 
similar: the fixed length constraint is best realized when 
more than a single irrep has condensed. So, for NVO 
and TMO as the temperature is lowered one encounters 
a second phase transition in which a second irrep appears. 
Within a low-order Landau expansion this phenomenon 
is described by a free energy of the form& 



F = \{T-T y )a 



-(T-T<)a| +u > cr 



+u < a\ + wa > a < , 



(48) 



where T> > T<. This system has been studied in de- 
tail by Bruce and Aharony^ For our purposes, the most 
important result is that for suitable values of the pa- 
rameters ordering in cr> occurs at T> and at some lower 
temperature order in <t< may occur. The application of 
this theory to the present situation is simple: we can (and 
usually do) have two magnetic phase transitions in which 
first one irrep and then at a lower temperature a second 
irrep condense. A question arises as to whether the con- 
densation of one irrep can induce the condensation of a 
second irrep. This is not possible because the two irreps 
have different symmetry. But could the presence of two 
irreps, T> and F< induce the appearance of a third irrep 
F3 at the temperature at which T< first appears? For 
that to happen would require that T™ (g> ® con- 
tain the unit representation for some values of n and m. 
This or any higher combination of representations is not 
allowed for the simple four irreps system like NVO. In 
more complex systems one might have to allow for such 
a phenomenon. 



III. APPLICATIONS 

In this section we apply the above formalism to a num- 
ber of multiferroics of current interest. 



A. Mn WO , 

MnW04 (MWO) crystallizes in the space group P2/c 
(#14 in Ref. l28f) whose general positions are given in 
Table [V] The two magnetic Mn ions per unit cell are at 
positions 

ri = (-,</,-) , T 8 = (i,l-y,|) . (49) 

The wavevector of incommensurate magnetic ordering 
isM^ q = (q x , 1/2,(7,) with q x a -0.21 and q z » 0.46) 



Er = (a;,y,z) 
Ir = (x, y, z) 



m v r = (x,y,z+ §) 
2 y r = (x,y,z + |) 



TABLE V: General Positions for space group P2/c 



and is left invariant by the identity and m y . We start by 
constructing the eigenvectors of the quadratic free energy 
(i. e. the inverse susceptibility matrix) . Here we use unit 
cell Fourier transforms to facilitate comparison with Ref. 
l39l Below X, Y, and Z denote integers (in units of lattice 
constants). When 



R/ + r / = (X,Y,Z) + n 



then 



(X + l,-Y-y,Zi 



(50) 



2 

(X,-Y 



1,Z-1)+t 2 . 
Then Eq. I|18|l gives the eigenvalue condition to be 



(51) 



27Tiq-[(2Y+l)j+fc] 



Um y )S a (ci,T 2 )e m+2 ^ 
AS Q (q,Ti) , 



(52) 



where i, x {fny) 
R / + 



<y( m v) = tz{m y ) 



-1. When 



Tf = (X,Y,Z)+T 2 

1 



(X 



2' 



Y 



i-y,z + -) 



(53) 



then 



R, 



{X+ 1 -,-Y-l-y,Z+\) 



(54) 



= (X,-Y-l,Z) + Tl , 
and Eq. IjlSJl gives the eigenvalue condition to be 

S'M,r 2 ) = Umy)S a (c l ,T 1 ) e 2 ^< 2Y +^ 

= £ a (m y )S a (q.,Ti)[-l] = A5 a (q, r 2 ) .(55) 

From Eqs. JS2J| and we get A = ±e I7r ^ z and 

S a (q.,T 2 ) = -[^ a (m y )/X\S a (q,n) . (56) 

So we get the results listed in Table IvTl 

So far the analysis is essentially the completely stan- 
dard one. Now we use the fact that the free energy is 
invariant under spatial inversion, even though that oper- 
ation does not conserve wavevector ^iSii We now deter- 
mine the effect of inversion on the n's. As will become 
apparent use of unit cell Fourier transforms makes this 
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Irrep 


Ti 


r 2 




Irrep 




r 2 


\(m y ) = 








\{m y ) = 


e inq z 






a ti x 


Q* Tlx 






* 

a t 


• * 

— 2<i r 


S(q,l) 


a*n y 


a*n y 




S(q,l) 


ia* s 


a*s 




a*n z 


a*n z 






a t 


—ia*t 




an x 


—an x 






ar 


iar 


S(q,2) 


—an y 


an y 




S(q2) 


—ias 


as 




an z 


—an z 






at 


iat 



TABLE VI: Allowed spin eigenfunctions for MWO (apart 
from an overall phase factor) before inversion symmetry is 
taken into account, where a — exp(— mq z /2). Here the n(q)'s 
are complex and we have taken the liberty to adjust the overall 
phase to give a symmetrical looking result. But these results 
are equivalent to Table II of Ref. 133 . 



analysis more complicated than if we had used actual 
position transforms. We use Eq. (|21|l to write 

XS(q, r = 1) = S(q, r = 2 )* e ~ 27 " e '^ +3 ' + * ; ) 

= bS(q,2)* , (57) 

where b = — exp[— 2iri(q x + q z )]. For T2 we get 

T[n x ,n v ,n z ] = [-n x ,n y ,-n z ]*b , (58) 

which we can write as 

ln a = b^ a (m v )n* a . (59) 

Now the free energy is quadratic in the Fourier spin coef- 
ficients, which are linearly related to the n's. So the free 
energy can be written as 



Fn 



n'Gn 



(60) 



where n = (n x , n y , n z ) is a column vector and G is a 3 x 3 
matrix which we write as 



G 



A a (3 
a* B 7 
13* 7* C 



(61) 



where, for Hermiticity the Roman letters are real and 
the Greek ones complex. Now we use the fact that also 
we must have invariance with respect to inversion, which 
after all is a crystal symmetry. Thus 



F 2 = [Zn]tG[In] • 
This can be written as 

F 2 = ^b€ a {m v )n a G a pb*a*t;j3(my)n* 



(62) 



a/3 



= )£c t {m v )n a G a p£p{my)n* . 

a/3 



(63) 



TABLE VII: As Table IVT1 except that here inversion sym- 
metry is taken into account. Here r, s, and t are real. All 
six components can be multiplied by an overall phase factor 
which we have not explicitly written. 



Thus we may write 



Fo 



= nt 



A -a (3 
—a* B —7 
13* -7* C 

A ~a* P* 
-a B -7* 
-7 C 



n , 



(64) 



where "tr" indicates transpose (so n tr is a row vector). 
Since the two expressions for F 2 , Eqs. l|50"|) and (jHU, 
must be equal we see that a — ia, (3 — b, and 7 = ic, 
where a. b, and c must be real. Thus G is of the form 



G 



A ia b 
—ia B ic 
b -ic C 



(65) 



where all the letters are real. This means that the critical 
eigenvector describing the long range order has to be of 
the form 



i Hz) 



= e l(t > 



(r,is,t) , 



(66) 



where r, s, and t are real. For T 2 we set e 1 ^ = —i. For Ti 
a similar calculation again yields Eq. I|66|l . but here we 
set — 1. (These choices are not essential. They just 
make the symmetry more obvious.) Thus we obtain the 
final results given in Table IVIII Lautenschlager et alSi 
say (just above Table II) "Depending on the choice of the 
amplitudes and phases ..." What we see here is that in- 
version symmetry fixes the phases without the possibility 
of a choice (just as it did for NVO). Note again that we 
have about half the variables to fix in a structure deter- 
mination when we take advantage of inversion invariance 
to fix the phase of the complex structure constants. 



1. Order Parameter 

Now we discuss the definition of the order parameter 
for this system. For this purpose we replace r by err, s 
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by crs. etc., with the normalization that 



r 2 + s 2 + t 2 = 1 



(67) 



Here the order parameter cr is complex because we al- 
ways have the freedom to multiply the wavefunction by 
a phase factor. (This phase factor might be "locked" by 
higher order terms in the free energy, but we do not con- 
sider that phenomenon here4£) We record the symmetry 
properties of the order parameter. With our choice of 
phases we have 



Zcr n (q) 
m y cr n (q) 
m y cr n (-q) 



= K(q)]* , 

= A(r„)er„(q) 

= A(r„)V„(-q) , 



(68) 



where <x n (q) is the complex- valued order parameter for 
ordering of irrep T n and A(r ra ) is the eigenvalue of m y 
given in Table IVm Now we write an explicit formula for 
the spin distribution in terms of the order parameters of 
the two irreps: 

S(R, t=1) = 2(7i (ni + t 1 k) cos(q ■ R + fa - nq z /2) 
+s h j sin(q • R + fa - irq z /2) 



-2a* 



{-r 2 i - t 2 k) sin(q • R + fa - nq z /2) 



+s 2 j cos(q • R + fa ~ , 



(69) 



S(R, t = 2) = 2(7i (ni + t 1 k) cos(q ■ R + fa + nq z /2) 
- sij sin(q • R + fa + nq z /2) 



-2a 2 



(r 2 i + t 2 k) sin(q • R + fa + nq z /2) 



+s 2 j cos(q • R + fa + irq z /2) 



One can explicitly verify that these expressions are con- 
sistent with Eq. (|68|> . Note that when only one of the 
order parameters (say a n ) is nonzero, we have inver- 
sion symmetry with respect to a redefined origin where 
fa = 0. For each irrep we have to specify three real 
parameters, ar n , <rs ni and at n and one overall phase fa 
rather than three complex-valued parameters had we not 
invoked inversion symmetry. 



B. TbMnOs 

Here we give the full details of the calculations for 
TbMnC>3 described in Ref. y. The presentation here 
differs cosmetically from that in Ref. 0. The space group 
of TbMn0 3 is Pbnm which is #62 in Ref. [H (although 
the positions are listed there for the Pnma setting). The 
space group operations for a general Wyckoff orbit is 
given m Table IVTTT1 In Table |IX| we list the positions 



Er = (x,y,z) 


2 x r=(x + ±,y + 




2z = (x,y,z+ |) 


2„ = (a; + i,y + i 




Jr = (x,y,z) 


m,r = (x+ \ ,y-V 




m z v = (x,y,z+ |) 


m y r= (x+ \,y+ i 


,z + \) 



TABLE VIII: General Positions for Pbnm. Notation as in 
Table □ 



Mn 


(1) = (0,|,0) 
(3) = (0,U) 


(2) = (1,0,0) 
(4) = (i0,i) 


Tb 


(5) = (x,y,i) 
(7) = (i,»,|) 


(6) = (x+i,y+|,|) 
(8) = (x + |,z/+i,i) 



TABLE IX: Positions of the Magnetic Ions in the Pbnm Struc- 
ture of TbMn0 3 , with x = 0.9836 and y = 0.08104i 



of the Mn and Tb ions within the unit cell and these are 
also shown in Fig. [3] 

To start we study the operations that leave invariant 
the wavevector of the incommensurate phase which first 
orders as the temperature is lowered. Experimentally 42 
this wavevector is found to be (0,q,0), with 34 q 



0.28(2tt/6). These relevant operators (see Table IvTTTf ) 
m x and m z . We follow the approach used for MWO, 
but use "actual location" Fourier transforms. We set 
R/ + Tt = r in order to use Eq. (|16|l and we need to 
evaluate 



A = exp^27rz(q/6) • [r - [mj 

>(2-K{qfb)j ■ [yj - [mj- 1 ^']) = ^ (71) 



exp 



(70) and 



A' = cxp^27ri(q/6) • [r - [m z ] : r 

= exp( 2m{q/b)3 ■ [yj - [m^yj] ) = 1 . (72) 



We list, in Table Ixl the transformation table of sublattice 
indices of TMO. 

Therefore the eigenvalue conditions for transformation 
by m x are 



S' a (<hTf) = Z a (m x )S a {c l ,T i )A=\(m x )S a (c l ,T f l73) 



and 



S'M^f) = £ a (rn z )S a (c l ,T i ) = X(m z )S a {q,T f ) (74) 

where £ x (m x ) = -£, y {m x ) = -£, z {m x ) = 1 and £ a (m z ) 
was defined in Eq. i|15|) . From these equations we see 
that X(m x ) assumes the values ±A and A(m z ) the values 
±1. Then solving the above equations leads to the results 
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z 




X 



FIG. 3: (Color online). Mn sites (smaller circles, on-line red) 
and Tb sites (larger circles, on-line blue) in the primitive unit 
cell of TbMnOs. The Tb sites are in the shaded planes at 
z — n ± | and the Mn sites are in planes z = n or = n + | , 
where n is an integer. The incommensurate wavevector is 
along the b axis. The mirror plane at z = 1/4 is indicated 
and the glide plane m x is indicated by the mirror plane at 
x = 3/4 followed by a translation (indicated by the arrow) of 
6/2 along the i/-axis. 



Ti 


T f (m x ) 


7-/(771*) 
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7 
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5 
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6 



TABLE X: Transformation table for sublattice indices of 
TMO under various operations. 

given in Tabic IXII (These results look different than 
those in Ref. because here the Fourier transforms are 
defined relative to the actual positions, whereas there 
they are defined relative to the origin of the unit cell.) 

Now, since the crystal is centrosymmetric, we take 
symmetry with respect to spatial inversion, I, into ac- 
count. As before, recall that X transports the spin to its 
spatially inverted position without changing the orienta- 
tion of the spin (a pseudovector). The change of position 
is equivalent to changing the sign of the wavevector in the 
Fourier transform and this is accomplished by complex 
conjugation. Since the Mn ions sit at centers of inversion 
symmetry, one has, for the Mn sublattices, 

XS(q,n) = S(q,n)* , (75) 
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TABLE XI: Spin functions (i. e. actual position Fourier co- 
efficients) within the unit cell of TMO for wavevector (0, q, 0) 
which are eigenvectors of m x and m z with the eigenvalues 
listed, with A = exp(— inq). All the parameters are complex- 
valued. The irreducible representation (irrcp) is labeled as in 
Ref. 0- Inversion symmetry is not yet taken into account. 
Note that the two Tb orbits have independent complex am- 
plitudes. 

where the second argument specifies the sublattice, as 
in Table IIXI In order to discuss the symmetry of the 
coordinates we define x\ = n a M , x-i = n b M , X3 = n c M and 
for irreps Ti and T^, and whereas 

for irreps T2 and T4, X4 = n^ 1 , X5 = n^ 2 , xq = n^, and 
Thus Eq. I|75(l gives 

2x n = a;* , n = 1, 2, 3 . (76) 

For the Tb ions X interchanges sublattices #5 and #7 
and interchanges sublattices #6 and #8. So we have 

XS(q,5) - S(q,7)* 

JS(q,6) - S(q,8)* . (77) 

Therefore we have 

2x4 = £5 , Zxq = . (78) 

Now we use the invariance of the free energy under X 
to write 

X,a;Y,{3 
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* nra^m 



^ ^ f^Jyi] G n7TL [Xx m ] , 



(79) 



where the matrix G is Hcrmitian. 

For irreps Ti and T3 the matrix G in Eq. (|79(l couples 
five variables. Equation 1|76[1 implies that the upper left 
3x3 submatrix of G is real. Equations (|7uT) and l(7S|) 

imply that G n ,4 — G§^ n and for n = 1,2, 3. We find that 
G assumes the form 



G = 



a b c a a* 

b d e /3 P* 

c e / 7 7* 

a* p* 7* 9 6 

a P 7 6* g 



(80) 



where the Roman letters are real valued and the Greek 
are complex valued. As shown in the appendix, the form 
of this matrix ensures that the critical eigenvector can be 
taken to be of the form 

ip = (nl I ,n b M ,n c M , n^n^*) = (r, s,t; p, p*) , (81) 

where the Roman letters are real and the Greek ones 
complex. Of course, because the vector can be complex, 
we should include an overall phase factor (which amounts 
to arbitrarily placing the origin of the incommensurate 
structure), so that more generally 



(82) 



For irreps Y2 and T4 the matrix G in Eq. (|79|) couples 
the seven variables listed just below Eq. (|75|l . Equations 



(|?6*|) and f?%|l imply that G n ^ = G*5, n and G n fi 
for n = 1, 2, 3. Therefore G assumes the form 



= Gt n 
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(83) 



where Roman letters are real and Greek are complex. 
As shown in the appendix, this form ensures that the 
eigenvectors are of the form 
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L Mi 11 Mi 
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l T1 , Iirp2 



(84) 



These results are summarized in Table IXIII Note that 
the use of inversion symmetry fixes most of the phases 
and relates the amplitudes of the two Tb orbits, thereby 
eliminating almost half the fitting parameters^ 
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TABLE XII: As Table 001 Apart from an overall phase <j>r 
for each irrep, inversion symmetry restricts all the manganese 
Fourier coefficients to be real and all the Tb coefficients to 
have the indicated phase relations. 



Order Parameters 



We now introduce order parameters cr n (q) = a n e 



for irrep T n in terms of which we can write the spin dis- 
tribution. For instance under T3 one has 
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S z {v,T2) 


= 2pa 3 cos{qy + 


h - <t> P ) 



(85) 



where we set p 
ized by 



p e *0p an( ^ ^- ne parameters are normal- 



+ s 2 + t 2 + \p\ 2 = 1 . 



(86) 



In Eq. (|85|l r = (x, y, z) is the actual position of the spin 
in question. From Table IXll one can obtain the symmetry 
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Er = 


(.x,y,z) 


2 x r = 


(x + \,y + i,z) 


2 Z = 


(x, y,z) 


2 y = 


{x + ±,z) 


Ir = 


(x,y, z) 


m x r = 


(x+±,y+ \,z) 


m z r = 


= (x,y,z) 


m y r = 


(x + |,f + \,z) 



TABLE XIII: As Table IXTTTI General Positions for Pbam. 



properties of the order parameters for each irrep. 
instance 



For 



m x <T 2 (q) 
m x cr 3 (q) 
m x (T 4 (q) 



and 



-Ao-i(q) 
-Aer 2 (q) 
-A<T 3 (q) 
-Acr 4 (q) 



I<r„(q) 



%ffi(q) 
m z o- 2 (q) 
m z cr 3 (q) 
m z <T 4 (q) 



;(q) 



+o-i(q) 

-o- 2 (q) 
CT 3 (q) 

-«r 4 (q)(87) 



Note that in contrast to the case of NVO, inversion 
symmetry does not fix all the phases. However, it 
again drastically reduces the number of possible mag- 
netic structure parameters which have to be determined. 
In particular, it is only by using inversion that one finds 
that the magnitudes of the Fourier coefficients of the two 
distinct Tb sites have to be the same. Note that if we 
choose the origin so that = (which amounts to re- 
naming the origin so that that becomes true), then we 
recover inversion symmetry (taking account that inver- 
sion interchanges terbium sublattice #3 and #1). One 
can determine that the spin structure is inversion invari- 
ant when one condenses a single representation. 

The result of Table 5 applies other manganates pro- 
vided their wavevector is also of the form (0, q y , 0). This 



includes YMn0.4i and HoMnO-, 



45.46 



Both these sys- 



tems order into an incommensurate structure at about 
T c w 42K. The Y compound has a second lower- 
temperature incommensurate phase, whereas the Ho 
compound has a lower-temperature commensurate phase. 



C. TbMn 2 5 

The space group of TbMn 2 5 (TM025) is Pbam (#55 
in Ref. l28f) and its general positions are listed in Table 
IXIHl The positions of the magnetic ions are given in 
Table IXTVl and are shown in Fig. 0] 

We will address the situation just below the ordering 
temperature of 43KiA2, We take the ordering wavevector 
to beiS to be (1,0, (?) withJi q m 0.306. (This may be 
an approximate value. 49 ) Initially we assume that the 
possible spin configurations consistent with a continu- 
ous transition at such a wavevector are eigenvectors of 
the operators m x and m y which leave the wavevector in- 
variant. We proceed as for TMO. We use the unit cell 




FIG. 4: (Color online). Two representations of TbM^Os. 
Top: Mn sites (on-line red) with smaller circles Mn 3+ and 
larger circles 4+ and Tb sites (squares, on-line blue) in the 
primitive unit cell of TbM^Os. The Mn +4 sites are in the 
shaded planes at z — n ± 8 with 8 w 0.25 and the Mn +3 sites 
are in planes z = n, where n is an integer. The Tb ions are in 
the planes z — n + i . The glide plane m x is indicated by the 
mirror plane at x = 3/4 followed by a translation (indicated 
by the arrow) of b/2 along the y-axis and similarly for the 
glide plane m y . Bottom: Perspective view. Here the Mn 3+ 
are inside oxygen pyramids of small balls and the Mn 4+ are 
inside oxygen octahedra. 



Fourier transforms and write the eigenvector conditions 
for transformation by m x as 

S a (cL,T f )' = &,(m s )S a (q,T i )e fa i< r '- R «> 

= A x ,5 a (q,r / ) , (89) 

where Tj and are respectively the sublattice indices 
and unit cell locations before transformation and Tf and 
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Mn 3+ 


(l) = (x,y,0) 
(3)(<E+i,y+i,0) 


(2) = (af,y,o) 

(4) = (x+|,y+i,0) 


Mn 4+ 


(5) = (|,0,2) 
(7) = (|,O,0) 


(6) = (0,i,^) 
(8) = (0,|,I) 


RE 


(9) = (x,y,i) 
(ii) = (x + i,r+i,i) 


(io) = (x,y,|) 
(i2) = (x + |,y + |,|) 



TABLE XIV: Positions of the magnetic ions of TbMn 2 O s in 
the Pbam structure. Here x = 0.09, y = -0.15, z = 0.25, 47 
X = 0.14, and Y = 0.174&. All these values are taken from 
the isostructural compound HoMn2 0s. 



R/ are those after transformation. The eigenvalue equa- 
tion for transformation by m y is 



S a (q,T f )> 



ta(my)S a (ti,T^'f-^ 
X y S a (q,T f ) . 



(90) 



If one attempts to construct spin functions which are 
simultaneously eigenfunctions of m x and m y one finds 
that these equations yield no solution. While it is, of 
course, true that the operations m x and m y take an 
eigenfunction into an eigenfunction, it is only for irreps 
of dimension one that the initial and final eigenfunc- 
tions are the same, as we have assumed. The present 
case, when the wavevector is at the edge of the Bril- 
louin zone is analogous to the phenomenon of "sticking" 
where, for nonsymmorphic space group (i. e. those hav- 
ing a screw axis or a glide plane) the energy bands (or 
phonon spectra) have an almost mysterious degeneracy 
at the zone boundary^ and the only active irrep has di- 
mension two. This means that the symmetry operations 
induce transformations within the subspace of pairs of 
eigenfunctions. We now determine such pairs of eigen- 
functions by a straightforward approach which do not 
require any knowledge of group theory. Here we explic- 
itly consider the symmetries of the matrix x 1 f° r the 
quadratic terms in the free energy which here is a 36 x 36 
dimensional matrix, which we write as 



yiixx) y^{xy) y^{xz) 

M txj/)t M to) M^^ 
M.( zx ^ M^ )t M^^ 



(91) 



where M^ afc '' is a 12 dimensional submatrix which de- 
scribes coupling between a-component and ^-component 
spins and is indexed by sublattice indices r and r' The 
symmetries we invoke are operations of the screw axes, 
m x and m y which conserve wavevector (to within a re- 
ciprocal lattice vector), and X, whose effect is usually 
ignored. To guide the reader through the ensuing cal- 
culation we summarize the main steps. We first ana- 
lyze separately the sectors involving the x, y, and z spin 
components. We develop a unitary transformation which 
takes M' Q ") into a matrix all of whose elements are real. 
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TABLE XV: Transformation table for sublattice indices with 
associated factors for TM025 under various operations, as 
defined by Eq. I|19jl . For m,, one has exp[iq ■ (R/ — Ri)] = 
1 for all cases and for m x m y I the analogous factor is +1 
in all cases and this operator relates S a (q, t) and S a (q, r)*. 
NOTE: This table does not include the factor of £ a (O) which 
may be associated with an operation. 

a) 4> = q ■ (R/ — Ri), as required by Eq. 1181 . 

b) ((>' = q • ("7% + Tf), as required by Eq. I2H . 



This fixes the phases within the 12 dimensional space 
of the a spin components within the unit cell (assuming 
these relations are not invalidated by the form of M( Q ^, 
with a 7^ P). The relative phases between different spin 
components is fixed by showing that the unitary transfor- 
mation introduced above leads to M> ay ) having all real- 
valued matrix elements and M^ 2 ) and M( yz ' having all 
purely imaginary matrix elements. The conclusion, then, 
is that the phases in the sectors of x and y components 
are coupled in phase and the sector of z components are 
out of phase with the x and y components. 



1. x Components 



As a preliminary, in Table IXVI we list the effect of 
the symmetry operations on the sublattice index. When 
these symmetries are used, one finds the 12 x 12 sub- 
matrix of M^ xx ' which couples only the x-components of 
spins assumes the form 
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(92) 
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where Roman letters are real quantities and Greek ones 
complex. (In this matrix the lines are used to separate 
different Wyckoff orbits.) The numbering of the rows and 
columns follows from Table IXTVl I give a few examples 
of how symmetry is used to get this form. Consider the 
term Ti, where 

Ti = xri^(-q,i)^(q,5) . (93) 

Using Table IXVl we transform this by m x into 

T[ = XlfiS*(-<L3)S x (sh6) , (94) 

which says that the 1,5 matrix element is equal to the 
3,6 matrix element. (Note that in writing down T[ we 
did not need to worry about £ a , since this factor comes 
in squared as unity.) Likewise if we transform by m y we 
get 



T[ = xr 6 1 [-5 a (-q s 4)][5 ai (q s 6)] 



(95) 



which says that the 1,5 matrix element is equal to the 
negative of the 4,6 matrix element. If we transform by 
m x m y we get 



T[ = xri[S*(-q,2)][-< ?x (q,5)] , 



(96) 



which says that the 1,5 matrix element is equal to the 
negative of the 2,5 matrix element. To illustrate the effect 
of X on T\ we write 



T[ = xrj[^(q,2)][-^(-q,7)] , 



(97) 



so that the 1,5 element is the negative of the 7,2 element. 
From the form of the matrix in Eq. (|92|) (or equivalcntly 
referring to Table IX*Xl in Appendix B), we see that we 
bring this matrix into block diagonal form by introducing 
the wavefunctions for S x (q, r) , 
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(98) 



The superscripts a, n on O label, respectively, the Carte- 
sian component and the column of the irrep accord- 
ing to which the wavefunction transforms. The sub- 
scripts m, t label, respectively, the index number of the 
wavefunction and the sublattice label. In this subspace 
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(99) 



where the coefficients are separated into real and imag- 
inary parts as V2a = a' + ia", \f2(3 = /3' + i(3" 
V27 = 7' + 17", and \/25 = 5' + iS" '. There are no 
nonzero matrix elements between wavefunctions which 
transform according to different columns of the irrep. 
The partners of these functions can be found from 



(100) 



so that, using Table IXVl and including the factor £ Q , we 
get 
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Within this subspace the matrix (n\M^ xx ' \m) is the same 
as in Eq. I|99|) because 

(nlm^M^mylm) = (n\M {xx) \m) . (102) 

These functions transform as expected for a two di- 
mensional irrep, namely, 











m x 

























(103) 



We will refer to the transformed coordinates of Eqs. 
(|98J) and IjlOljl as "symmetry adapted coordinates." The 
fact that the model-specific matrix that couples them is 
real, means that the critical eigenvector is a linear com- 
bination of symmetry adapted coordinates with real co- 
efficients. 



2. y Components 

The 12 x 12 matrix M^ yv ^ coupling y components of 
spin has exactly the same form as that given in Eq. 19211 , 
although the values of the constants are unrelated. This 
is because here one has ^ = 1 in place of £ x = 1. There- 
fore the associated wavefunctions can be expressed just 
as in Eqs. I|98|l and (|1C)1|) except that all the superscripts 
are changed from x to y and t now labels S y (q, t). How- 
ever, the transformation of the y components rather than 
the x components, requires replacing £ x by £ y which in- 
duces sign changes, so that 
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(104) 
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We want to construct wavefunctions in this sector which 
transform just like the x components, so that they can 
be appropriately combined with the wavefunctions for 
the ^-components. In view of Eq. I|103|) we set 



n(a,i) = n(x,2) n (y,2) _ n ( x ,i) 

n.r n.r ' n.r 7i,r 



(105) 



So the coefficients for On are given by Eq. I|101|) and 

those for O^' 2 ' by Eq. (|98|) . These wavefunctions are 
constructed to transform exactly as those for the x com- 
ponents. 

3. z Components 

Similarly, we consider the effect of the transformations 
of the z components. In this case we take account of the 
factor £ z to get 



rn. 



" (*,1) " 
















" (*,a) 









(106) 



We now construct wavefunctions in this sector which 
transform just like the x components. In view of Eq. 
(|TU5|> we set 



rt(z,i) _ n (x,2) r>(z,2) _ _ n (x,i) 



So the coefficients for oi 2 ' 1 ^ are given by Eq. I|101|) 

and those for 0„ ' are the negatives of those of Eq. 
(|98|l . These wavefunctions are constructed to transform 
exactly as those for the x components. 



(107) 



4- The Total Wavefunction and Order Parameters 

Now we analyze the form of M( ab ) of Eq. lj!?Tl) for a^b, 
using inversion symmetry. To do this it is convenient to 
invoke invariance under the symmetry operation m x m y I 
whose effect is given in Table IXVl We write 



m x m y lS a (q,T) = € a (m x )£ 

xS a (q,ftr)* , (108) 

where TZr = t for r ^ 5, 6, 7, 8, otherwise TZt = t ±2 
within the remaining sector of t's and a (and later b) 
denotes one of x, y, and z. Thus 

T = 5 (q, TfM^) S 6 (q,r') 

= [m x m y lS a (q, r)]*Mffi [m x m y lS b (q, t')] 

= C ab S a (q,TlT)M^ S b (<i,KT'y , (109) 

where 

Cab = £a(rn x )£a(m y )£ b (m x )£, (m y ) . (110) 



From this we deduce that 

r (ba 



(ab) 



or, since M is Hermitian that 



M^) = C ab \M^ b \ „ , , 

tt' uu 7c l t.1Z 1 r / 



(111) 



(112) 



Now we consider the matrices M( afc -', in the symmetry 
adapted representation where 

M (ab) = X^\o a P]*M (a> ! ) O bp , 

±v± n,m / j V^nrl tt' tut' 

= J2 C ^° 

= c ab Y,[o 



-\ap "1 * 
nr J 



o bp , 



J 7iUt\ 



M (ab) 

T.T 



O h *nr>- (H3) 



There are no matrix elements connecting p and p' ^ p 
and the result is independent of p. One can verify from 
Eqs. lEHl and (|TUT|) that 



O ap 



[o 



ap 1 ^ 
n.T\ i 



(114) 



so that 



= C o6 [M^]* . 



We have that C xy = —C xz 



(115) 
1, so that all the 



elements of M^*) are real and all the elements of ~My xz > 
and are imaginary. Thus apart from an over all 

phase for the eigenfunction of each column, the phases of 
all the Fourier coefficients are fixed. What this means is 
that the critical eigenvector can be written as 

2 6 / 

i' = E °v E rnxOfrri + r ny O^ 

p—1 n—1 \ 



(116) 



where the r's are all real-valued and are normalized by 

6 



EE^ 2 = !> 



(117) 



n— 1 a 



and cr p are arbitrary complex numbers. Thus we have 
the result of Table IXVll 
The order parameters are 



a x = are* , a 2 = a 2 e^ 



(118) 



Neither the relative magnitudes of &\ and <x 2 nor their 
phases are fixed by the quadratic terms within the Lan- 
dau expansion. Note that the structure parameters of Ta- 
ble IXVII are determined by the microscopic interactions 
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TABLE XVI: Normalized spin functions (i. e. Fourier co- 
efficients) within the unit cell of TbM^Os for wavevector 
(5,0,9). Here z a — (rz a + irA a )/j2. All the r's are real 
variables. The wavefunction listed under ai {02) transforms 
according to the first (second) column of the irrep. The ac- 
tual spin structure is a linear combination of the two columns 
with arbitrary complex coefficients. 



which determine the matrix elements in the quadratic 
free energy. (Since these are usually not well known, one 
has recourse to a symmetry analysis.) The direction in 
<xi-er 2 space which the system assumes, is determined by 
fourth or higher-order terms in the Landau expansion. 
Since not much is known about these terms, this direc- 
tion is reaosonably treated as a parameter to be extracted 
from the experimental data. We use Table IXVll to write 
the most general spin functions consistent with crystal 
symmetry as 

S(R, 1) = err (ri x i + rijj) cos(q • R + fa) 
+ri z fcsin(q ■ R + fa) 



-02 



(r 2x i + r 2y j] cos(q • R + fa) 



+r 2z fcsin(q ■ R + fa) 
S(R, 2) = eri {r 2x i + r 2y j) cos(q • R + fa ) 
-r 2z fcsin(q ■ R + fa) 
+<J2 {ri x i + ri y j]cos{q-R + fa) 



S(R 6 ) 



S(R, 3) = eri (n x i-r ly j)cos(q-R+fa) 
— r l2 /csin(q • R + cf)i) 



"02 



(~r 2x i + r 2y j] cos(q • R + fa) 



-fr 2z fcsin(q ■ R+ cj) 2 ) 



S(R,4) = a x 



(r 2x i - r 2y j) cos(q • R + fa) 
r 2z k s\n{q-R + fa) 



-02 



{-r\ x i + n y j] cos(q • R + fa) 



-r iz A:sin(q • R + fa) 



S(R,5) = tj 1 



(z' x i — z'j — z" z k) cos(q • R + fa) 



+{z"i - z"j + z' z k) sin(q ■ R + fa] 



-02 



(—z' x i + z' j — z" z k\ cos(q • R + fa) 



+ (—z"i + z''j + z' z k) sin(q • R - 



01 



{z'j, + z. ' j + z' z k) cos(q • R + fa) 



+(z"i + z'yji — z' z k) sin(q • R + fa) 



-o 2 



(z' x i + z. j — z'Jk] cos(q • R + fa) 



+ {z'J + zfj + z' z k) sin(q ■ R + <t>; 
S(R, 7) = o 1 {z'j - z' y j + z' z k)cos{q-R + fa) 
+ {-z"i + z"j + z' z k) sin(q • R + fa) 



+cr 2 {-z'j + z y j + z"k] cos(q • R + fa) 
+ {z"i - z'yj + z' z k) sin(q ■ R + fa) 
S(R 8 ) = <7\ {z'J + z'yj - z' z r k)cos{q-R + fa) 
+ {-zfi - zfj - z' z k) sin(q • R + fa) 



-02 



{z'j + Zyj + z" z k\ cos(q • R + fa) 



+(-z'fi - zf + z' z k) sin(q • R + fa) 



S(R,9) = a x 



{r 5x i + r 5y j) cos(q • R + fa) 
r 5z fcsin(q- R + fa) 



+cr 2 



[r 6x i + r 6y j] cos(q -R+fa) 



S(R, 10) = o-i 



-r 6z £;sin(q -R+fa) 

{r 6x i + r 6y j) cos(q -R + fa) 
r 6z &sin(q -R + fa) 



-r lz fcsin(q -R + fa) 
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+c 2 (r 5x i + r 5y j]cos(q-'R + (l)2) 

-r 52 fcsin(q- R + fa) 
S(R, 11) = o\ (r 5x i-r 5y j)cos{c[-R + fa) 

-r 5z fcsin(q-R + </)i) 

+ct 2 (-r 6x i + r 6y j] cos(q • R + fa) 

+r 6z fcsin(q • R + fa) 
S(R, 12) = <Ti (r 6x i-r 6y j)cos{c[-R + fa) 

+r 6z fcsin(q ■ R+ fa) 

+cr 2 (-r 5x i + r 5y j] cos(q • R + fa) 

-r 5z k sin(q-R + fa) 



(119) 



In Table IXVll the position of each spin is R + t„, where 
the t are listed in Table IXTVl and R is a Bravais lattice 
vector. The symmetry properties of the order parameters 
are 
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(120) 



We now check a few representative cases of the above 
transformation. If we apply m x to S^q, 1) we do not 
change the signs of the x component but do change the 
signs of the y and z components. As a result we get 
5(q, 3) except that a y has changed sign, in agreement 
with the first line of Eq. (|120fl . If we apply m y to 5(q, 1) 
we do not change the sign of the y component but do 
change the signs of the x and z components. As a result 
we get 5(q, 4) except that now cr\ is replaced by <T2 and 
cr 2 is replaced by «Ti, in agreement with the second line 
of Eq. (|120fl . When inversion is applied to 5(q, 1) we 
change the sign of R but not the orientation of the spins 
which are pseudovectors. We then obtain <S(q, 2) provid- 
ing we replace <Xi by erj and er 2 by cr\, in agreement with 
the last line of Eq. ((HUj) . 



5. Comparison to Group Theory 

Here I briefly compare the above calculation to the one 
using the simplest formulation of representation theory. 
The first step in the standard formulation is to find the 
irreps of the group of the wavevector. The easiest way 
to do this is to introduce a double group having eight 
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m y = 






m x m y 









-1 


-1 
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elements (see Appendix B) since we need to take account 
of the operator m y = — E. (This is done in Appendix B.) 
From this one finds that each Wyckoff orbit and each spin 
component can be considered separately (since they do 
not transform into one another under the operations we 
consider). Then, in every case the only irrep that appears 
is the two dimensional one for which we set 
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Indeed, one can verify that the functions in the sec- 
ond (third) column of Table IXVII comprise a basis vec- 
tor for column one (two) of this two dimensional irrep. 
One might ask: "Why have we undertaken the ugly de- 
tailed consideration of the matrix for i 7 ^?" The point is 
that within standard representation theory all the vari- 
ables in Table IXVT1 would be independently assigned ar- 
bitrary phases. In addition, the amplitudes for the Tb or- 
bits (sublattices #5, #6 and sublattices #7, #8) would 
have independent amplitudes. To get the results actually 
shown in Table IXVT1 one would have to do the equivalent 
of analyzing the effect of inversion invariance of the free 
energy. This task would be a very technical exercise in 
the arcane aspects of group theory which here we avoid 
by an exercise in algebra, which though messy, is ba- 
sically high school math. I also warn the reader that 
canned programs to perform the standard representation 
analysis can not always be relied upon to be correct. It is 
worth noting that published papers dealing with TM025 
have not invoked inversion symmetry. For instance in 
Ref. ^3 one sees the statement "As in the incommensu- 
rate case [3], each of the magnetic atoms in the unit cell is 
allowed to have an independent SDW, i. e., its own am- 
plitude and phase," and later on in Ref. |4f| "all phases 
were subsequently fixed ... to be rational fractions of 7r." 
Use of the present theory would eliminate most of the 
phases and would relate the two distinct Mn 4+ Wyckoff 
orbits (just as happened for TMO). 

Finally, to see the effect of inversion on a concrete level 
I analyze the situation within the Mn 3+ orbit and con- 
sider only the x components of spin. The inverse sus- 
ceptibility would then be the upper left 4x4 submatrix 
shown in Eq. i|92|) . Had we not used inversion symmetry 
this submatrix would be the same except that one would 
have had 
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IT 



(122) 



where r is real. One can verify that when r = the 
eigenvectors can be taken to have only real components, 
whereas when r / 0, the eigenvectors are complex with 
relative phases dependent on the value of r, 



6. Comparison to YMn^ Os 

YMn 2 5 (YM025) is isostructural to TM025, so its 
magnetic structure is relevant to the present discussion. 
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I will consider the highest temperature magnetically or- 
dered phase, which appears between about 20K and 45K. 
In this compound Y is nonmagnetic and in the higher- 
temperature ordered phase q z = 1/4, so the system is 
commensurate. But since the value of q z is not special, 
the symmetry of this state is essentially the same as that 
of TM025. Throughout this subsection the structural 
information is taken from Fig. 2 of Ref. [nj. (The up- 
permost panel is mislabeled and is obviously the one we 
want for the highest temperature ordered phase.) 





> 4 


> * 







FIG. 5: (Color online). Top: The spin structure of the Mn 3+ 
ions in YM^Os (limited to one a-b pane), taken from Fig. 
2 of Ref. ITU The sublattices are labeled in our convention. 
Bottom left: the spin structure of after transformation by m x . 
Bottom right: spin structure of TbMn2 0s after transforma- 
tion by m y . 



From Fig. [S] we see that the spin wavefunction is an 
eigenvector of m x with eigenvalue —1. So this structure 
must be that of the second column of the irrep. In accor- 
dance with this identification one sees that the initial 
wavefunction is orthogonal to the wavefunction trans- 
formed by m y (since this transformation will produce a 
wavefunction associated with the first column). Refer- 
ring to Eq. (|119fl . one sees that to describe the pattern 
of Mn 3+ spins one chooses 

<7i=0, r 2x = -ri x « 0.95 , 

n y = -r 2y w 0.3 . (123) 

The point we make here is that <j\ =0. Although the 
values of these order parameters were not given in Ref. 
ITU it seems clear that in the lower temperature phase 
the order parameters must be comparable in magnitude. 



D. CuFe0 2 

The space group of CuFe0 2 (CFO) is 53 R3m (#166 
in Ref. l28i) and its general positions are given in Table 

IxvTTl 

We are interested in structures that can appear for 
general wavevectors of the type q = (q, q, 0) (in crystallo- 
graphic notation) , in other words for wavevectors parallel 



Er = (x,y,z) 
m 3 r = (y,x,z) 


3r = (z,x,y) 
m 2 r = (z,y,x) 


3 2 r = (y, z, x) 
mir = (x, z, y) 


It = (x, y, z) 
lm 3 r = {y,x,z) 


J3r = (z, x, y) 
Im 2 r = (z,y,x) 


T3 2 r = (y,z,x) 
Imir = (x,z,y) 



TABLE XVII: General Positions for R3m, with respect to 
rhombohedral axes. Here "3" denotes a three- fold rotation 
and m„ labels the three mirror planes which contain the three- 
fold axis. 




q 

X 



FIG. 6: The lattice of magnetic Fe ions in CFO. Here I 
show sections of three adjacent triangular lattice layers. The 
wavevector lies along the x-axis. The dashed line indicates 
that the central site lies directly above (below) the center of 
gravity, indicated by a dot, of a triangle of the layer below 
(above) it. 



to a nearest neighbor vector of the triangular plane of Fe 
ions. The only operation (other than the identity) that 
conserves wavevector is 2 X a two-fold rotation about the 
axis of the wavevector. Clearly, the Fourier component 
m z (q) obeys 

2 x m x (q) = m x (q) (124) 
and we call this irrep #1. For irrep #2 we have 



2 x m y {q) 
2 x m z (q) 



-rrij,(q) 
-m z (q) 



(125) 



As before inversion fixes the phase of these coefficients to 
be the same (within a given irrep), so one has 



Ol 



and 



rr>„ 



a 2 r , m z 



er 2 s 



(126) 



(127) 
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where r 2 + s 2 = 1 and er n (±|q|) = ^e*. Thus, if 
both irreps are present, we would have 

m x (r) — 2ai cos(qx + 

m y (r) = 2a 2 r cos(qx + <p 2 ) 

m z (r) = 2&2S cos(qx + $2) ■ (128) 

We have the transformation properties 

2 x cri — <T\ , 2 x cr 2 = —cr 2 
Jcr 1 = [ax]* , Xct 2 = [<r 2 ]* . (129) 

For future reference it is interesting to note that at zero 
applied electric and magnetic fields the free energy must 
be invariant under taking either <j\ or o 2 into its nega- 
tive. To see this write the free energy as an expansion in 
powers of the order parameters: 

F = ^ afci m „<Ti<Tj i <T™<T2™ . (130) 

To be time reversal invariant, the total number of powers 
(k + 1 + m + n) must be even. Also, to be invariant under 
2 X the total number of powers of u 2 (m + n) has to be 
even. So both k + 1 and m + n are even. But wavevector 
conservation implies that I + n = k + m (assuming we 
have a truly incommensurate phase). Thus I = k + 2a 
and n = m — 2a. Then 

F = J2 a kimn^ k+2a r a 2 2 m - 2 ° r e 2 M^-0i) . ( 131 ) 

E. Discussion 

1. Summary of Results 

In Table IXVIIII we collect the results for various mul- 
tiferroics. 



2. Effect of Quartic Terms 

As we now discuss, the quartic terms in the Landau ex- 
pansion can have significant qualitative effects^ In gen- 
eral, the quartic terms are the lowest order ones which 
favor the fixed length spin constraint, a constraint which 
is known to be dominant at low temperature^ How this 
constraint comes into play depends on what state is se- 
lected by the quadratic terms. For instance, in the sim- 
plest scenario when one has a ferromagnet or an antiferro- 
magnet, the instability is such (see Fig. that ordering 
with uniform spin length takes place. Thus, as the tem- 
perature is lowered within the ordered phase, the order- 
ing of wavevectors near q — for the ferromagnet (near 
q = 7r for the antiferromagnet) which would have be- 
come unstable if only the quadratic terms were relevant, 
is strongly disfavored by the quartic terms. In the sys- 
tems considered here the situation is quite different. For 



instance, in NVO^ TMOr and MWOS the quadratic 
terms select an incommensurate structure in which the 
spins are aligned along an easy axis and their magni- 
tudes are sinusoidally modulated. As the temperature is 
lowered the quartic terms lead to an instability in which 
transverse spin component break the symmetry of the 
longitudinal incommensurate phase. This scenario ex- 
plains why the highest-temperature incommensurate lon- 
gitudinal phase becomes unstable to a lower-temperature 
incommensurate phase which has both longitudinal and 
transverse components which more nearly conserve spin 
length. 

To see this result formally for NVO, TMO, or MWO, 
let <r> (<x<) be the complex valued order parameter for 
the higher-temperature longitudinal (lower-temperature 
transverse) ordering. The fourth order terms then lead 
to the free energy as 

F = a (T-r > )|<T > | 2 + 6(T-T < )| < T < | 2 
+A(|<x > | 2 + |,x < | 2 ) 2 + J B|<x > <x < | 2 
+C[(<x < <x^) 2 + «<x > ) 2 ] , (132) 

where A, B, and C are real. That C is real is a result 
of inversion symmetry, which, for these systems leads to 
Tcr n = <r*. The high-temperature representation does 
allow transverse components and could, in principle, sat- 
isfy the fixed length constraint. In the usual situation, 
however, the exchange couplings are nearly isotropic and 
this state is not energetically favored. If the higher tem- 
perature structure is longitudinal, then B will surely be 
negative, whereas if the higher temperature structure 
conserves spin length B will probably be positive. By 
properly choosing the relative phases of the two order 
parameters the term in C always favors having two ir- 
reps. So the usual scenario in which the longitudinal 
phase becomes unstable relative to transverse ordering is 
explained (in this phenomenology) by having B be neg- 
ative, so that the discussion after Eq. (|48|l applies. 

To finish the argument it remains to consider the term 
in C, which can be written as 

5F 4 = 2Ca 2 < a 2 > cos(20< - 2</»>) , (133) 

where again we expressed the order parameters as in Eq. 
(|43|1 . Normally, if two irreps are favored, it is because 
together they better satisfy the fixed length constraint. 
What that means is that when spins have substantial 
length in one irrep, the contribution to their spin length 
from the second irrep is small. In other words, the two 
irreps are out of phase and we therefore expect that to 
minimize SF4 we do not set </>< = 4>>, but rather 

0<=0>±7r/2. (134) 

In other words, we expect C in Eq. I|133|) to be positive. 
The same reasoning indicates that the fourth order terms 
will favor <j> 2 - <j>i = tt/2 in Eq. JTHJl for CFO. 

For all of these systems which have two consecutive 
continuous transitions one has a family of broken sym- 
metry states. At the highest temperature transition one 
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Phase 


T<(K) 


T>(K) 


q 


Irreps 


Refs. 


FE? 


Refs. 


NVO (HTI) 


6.3 


9.1 


(q,0,0) 


r 4 


6.33 


No 


4.6 


NVO (LTI) 


3.9 


6.3 


(q,0,0) 


r 4 + ri 


6,33 


1 b 


4,6 


TMO (HTI) 


28 


41 


(0,9,0) 


r 3 


3.42 


No 


^ 


TMO (LTI) 




28 


(0,?,0) 


r 3 + r 2 


3 


II c 


2 


TbMn 2 5 (HTI) 


38 


43 




f(b) 


10.49 


No 


58 


TbMn 2 Os (LTI) 


33 


38 


(1,0, q) 


p(c) 


10,49 


1 & 




YM112O5 (C) (d) 


23 


45 


(|,o, i) 


p(b) 


11 


II b 


& 


YMn 2 5 (IC) 




23 


(« 5,0,9) 




11 


II & 


58 


CFO (c) (HTI) 


10 


14 


(9,9,0) 


Ti 


£1 


No 




CFO (LTI) 


0? 


10 


(5,9,0) 


ri +r 2 


off) 


_L c 


9 


MWO 


12.7 


13.2 


(9*, |,g z ) 


r 2 


39 


No 


40 


MWO 


7.6 


12.7 


(9k, §,9*) 


r 2 + ri 


39 


II & 


40 



TABLE XVIII: Incommensurate Phases of various multiferroics. Except for CFO each phase is stable for zero applied magnetic 
field for T< < T < T> . When T< = it means that the phase is stable down to the lowest temperature investigated. We give 
the incommensurate wavevector and the associated irreducible representations in the notation of our tables. In the column 
labeled "FE?" if the system is ferroelectric we give the direction of the spontaneous polarization, otherwise the entry is "No." 

a) At the highest temperature the value of q x might not be exactly 1/2. 
b) The irrep is the two dimensional one (see Appendix B). In the HTI phase only one basis vector is active, 
c) The irrep is the two dimensional one (see Appendix B). In the LTI phase both basis vectors are active. 

d) This phase commensurate, 
e) Data for CuFe0 2 is for H ^ 8T. 
f ) The magnetic structure was inferred from the existence of ferroelectricity. 



has spontaneously broken symmetry which arbitrarily se- 
lects between er> and — cr>. (This is the simplest scenario 
when the wavevector is not truly incommensurate.) Inde- 
pendently of which sign is selected for the order parame- 
ter <r>, one similarly has a further spontaneous breaking 
of symmetry to obtain arbitrarily either ia < or —io < . 
(Here, as mentioned, we assume a relative phase 7r/2 for 
<r < . In this scenario, then, there are four equivalent low 
temperature phases corresponding to the choice of signs 
of the two order parameters. 

The cases of TM025 and YM025 are different from 
the above because they have two order parameters from 
the same two-dimensional irrep and which therefore are 
simultaneously critical. At quadratic order one has SU2 
symmetry, but this is broken by quartic terms in the 
free energy. The symmetry of these is such that in the 
representation of Eq. (|120fl one has 

F = a{T-T c )[\cr l \ 2 + \cT 2 \ 2 } 

+ A{\CT 1 \ 2 + \CT 2 \ 2 ) 2 +B\ ( T 1 CT 2 \ 2 

+C [{rr l( r*) 2 + Krr 2 ) 2 ] . (135) 

Terms odd in cr 2 are not allowed according to Eq. (|120fl . 
Also wavevector conservation indicates that two variables 
must be at wavevector q and two at wavwevector — q. 
Also A, B, and C are real. That C is real is a result of 
symmetry under m y , as in Eq. (|120fl . Here the fourth 
order anisotropy makes itself felt as soon as the ordered 
phase is entered. One can see that the phase difference 



between <Ti and cr 2 depends on the term proportional 
to C. Since the fixed spin length constraint favors this 
phase difference to be 7r/2, we intuit that C is positive, 
so that er 2 = ±icri. In any case, after minimizing with 
respect to the relative phase of the two order parameters, 
one gets the sum of the B and C terms as 

6F 4 = {B-\C\)\<r 1( r 2 \ 2 . (136) 

If B — |C| is positive, then either <Xi = or <x 2 = 0. 
In the former case the state is odd under m x and in the 
latter case even under m x . If B — |C| is negative, then 
the quartic terms favor 

ki| = ka|. (137) 

It is amusing that the quartic terms in the Landau ex- 
pansion dictate that these are the two allowable scenarios 
unless one admits to having a multicritical point where 
B-\C\ =0. 

Now we first consider YM025 in its higher temper- 
ature commensurate (HTC) ordered phase. For it addi- 
tional fourth order terms occur because 4q is a reciprocal 
lattice vector, but these are not important for the present 
discussion. Here the analysis of Ref. ^ indicates (see the 
discussion of our Fig. |SJ that only a single order parame- 
ter condenses in the HTC phase. This indicates that en- 
ergetics must favor positive B— \C\ in this case. The ques- 
tion is whether B — |C| is also positive for TM025. As 
we will see in the next section one has ferroelectricity un- 
less the magnitudes of the two order parameters are the 
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same. For YM025 the HTC phase is ferroelectric and the 
conclusion that only one order parameter is active com- 
ports with this. However, for TM025 the situation is not 
completely clear. Apparently there is a region such that 
one has magnetic ordering without ferroelectricity>i2i&& If 
this is so, then TM025 differs from YM025 in that its 
high temperature incommensurate phase has two equal 
magnitude order parameters. 

IV. MAGNETOELECTRIC COUPLING 

Ferroelectricity is induced in these incommensurate 
magnets by a coupling which is somewhat similar to that 
for the so-called "improper ferroelectrics."— To see how 
such a coupling arises within a phenomenological pic- 
ture, we imagine expanding the free energy in powers of 
the magnetic order parameters which we have studied in 
detail in the previous section and also the vector order 
parameter for ferroelectricity which is the spontaneous 
polarization P, which, of course, is a zero wavevector 
quantity. If we had noninteracting magnetic and elec- 
tric systems, then we would write the noninteracting free 
energy, F non as 

p — i Vv- 1 p 2 

a 

+ ^E a r(r-T r )|a r (q)| 2 +0(a 4 ) (138) 
r 

The first term describes a system which is not close to 
being unstable relative to developing a spontaneous po- 
larization (since in the systems we consider ferroelectric- 
ity is induced by magnetic ordering) . The magnetic terms 
describe the possibility of having one or more phase tran- 
sitions at which successively more magnetic order param- 
eters become nonzero. As we have mentioned, the sce- 
nario of having two phase transitions in incommensurate 
magnets is a very common onej^i and such a scenario 
is well documented for both NVO^ and TMOi2i A 
similar phenomenological description of second harmonic 
generation has invoked the necessity of having simulta- 
neously two irrepsS Below we will indicate the existence 
of a term linear in P, schematically of the form —XM 2 P, 
where A is a coupling constant about which not much be- 
yond its symmetry is known. One sees that when the free 
energy, including this term, is minimized with respect to 
P one obtains the equilibrium value of P as 

(P) = xe^M 2 . (139) 

A. Symmetry of Magnetoelectric Interaction 

We now consider the free energy of the combined mag- 
netic and electric degrees of freedom which we write as 

F = F Don + F int . (140) 



In view of time reversal invariance and wavevector con- 
servation, the lowest combination of M(q)'s that can ap- 
pear is proportional to M a (—q)Mp(q). So generically 
the term we focus on will be of the form 

F^t = Ec Q/37 M Q (q)Af (3 (-q)F 7 , (141) 

where a, /?, and 7 label Cartesian components. But, as 
we have seen in detail, the quantities M a (q) are linearly 
related to the order parameter or (q) , associated with the 
irrep T. Thus instead of Eq. (|141f) we write 

fint = Yl ^rr' 7 ^r(qkrv(q)*P 7 • (142) 

r,r',7 

The advantage of this writing the interaction in this form 
is that it is expressed in terms of quantities whose sym- 
metry is manifest. In particular, the order parameters 
we have introduced have well specified symmetries. For 
instance it is easy to see that for most of the systems stud- 
ied here, magnetism can not induce ferroelectricity when 
there is only a single representation present £4 This fol- 
lows from the fact that for NVO and TMO, for instance, 

ZK\ 2 = Wn\ 2 , (143) 

as is evident from Eq. (|47|l . The interpretation of this is 
simple: when one has one representation, it is essentially 
the same as having a single incommensurate wave. But 
such a single wave will have inversion symmetry (to as 
close a tolerance as we wish) with respect to some lat- 
tice point. This is enough to exclude ferroelectricity. So 
the canonical scenario^ is that ferroelectricity appears, 
not when the first incommensurate magnetic order pa- 
rameter condenses, but rather when a second such or- 
der parameter condenses. Unless the two waves have the 
same origin, their centers of inversion symmetry do not 
coincide and there is no inversion symmetry and hence 
ferroelectricity will occur. One might ask whether or not 
the two waves (i. e. two irreps) will be in phase. The 
effect, discussed above, of quartic terms is crucial here. 
The quartic terms typically favors the fixed length spin 
constraint. To approximately satisfy this constraint, one 
needs to superpose two waves which are out of phase. 
Indeed the formal result, obtained below, shows that the 
spontaneous polarization is proportional to the sine of 
the phase difference between the two irreps 1* We now 
consider the various systems in turn. 

B. NVO, TMO, and MWO 

We now analyze the canonical magneto-electric inter- 
action in the cases of NVO, TMO, ad MWO. These cases 
are all similar to one another and in each case the order 
parameters have been defined so as to obey Eq. 
This relation indicates that if we may choose the origin of 
the incommensurate system so that the phase of the order 
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parameter at the origin of a unit cell is arbitrarily close 
to zero. When this phase is zero, the spin distribution of 
this irrep has inversion symmetry relative to this origin. 
In the case when only a single irrep is active, this sym- 
metry then indicates that the magnetic structure can not 
induce a spontaneous polarization^ As mentioned, in the 
high temperature incommensurate phases of NVO, TMO, 
and MWO only one irrep is present, and this argument 
indicates that the magneto-electric interaction vanishes 
in agreement with the experimental observation^*^ that 
this phase is not ferroelectric. 

We now turn to the general case when one or more 
irreps are present 4*£*§*£ We write the magneto-electric 
interaction as 

= o H P -> [^rT' 7 <Tr(q)0- r ,(q)* 

7IT' 

+^rT 7 oT'(q)oT(q)*] • (144) 

For this to yield a real value of P we must have Hermitic- 
ity: Arr'-f = Ap T . In addition, because this is an ex- 
pansion relative to the state in which all order parameters 
are zero, this interaction has to be inversion under all op- 
erations which leave this "vacuum" state invariant ,22^ 
In other words this interaction has to be invariant un- 
der inversion (which takes P 7 into — P 7 ). This condition, 
together with Eq. (|47|l indicates that .Arr'7 = —Ap'r-/- 
This condition taken in conjunction with Hcrmiticity in- 
dicates that Arr'-y is pure imaginary. Thus 

F int = ~ ^ P 7 r IT'7 [0T(q)0T'(q)* 

7IT' 

-^r(q)VHq)] . (145) 

Since usually we have only two differenet irreps, we write 
this as 

Pnt = y^r 7 P 7 g>g< sin(<^> - <ft<) . (146) 
1 

where r 1 is real. The fact that the result vanishes when 
the two waves are in phase is clear because in that case 
one can find a common origin for both irreps about which 
one has inversion symmetry. In that special case one has 
inversion symmetry and no spontaneous polarization can 
be induced by magnetism. The above argument applies 
to all three systems, NVO^ TMO^ and MWO. As we 
will see in a moment, it is still possible for inversion sym- 
metry to be broken and yet induced ferroelectricity not 
be allowed. 

We can also deduce the direction of the spontaneous 
polarization by using the transformation properties of the 
order parameters, given in Eq. @S). We start by ana- 
lyzing the experimentally relevant cases at low or zero 
applied magnetic field. For NVO the magnetism in the 
lower temperature incommensurate phase is described^**^ 
by the two irreps T4 and V\. One sees from Eq. I|4t)|) 
that the product a\u^ is even under m z and odd un- 
der 2 X . For the interaction to be an invariant, P 7 



has to transform this way also. This implies that only 
the b-component of the spontaneous polarization can be 
nonzero, as observed^ For TMO the lower temperature 
incommensurate phase at low magnetic field is described*^ 
by irreps T3 and T2- From Table IXTfl we see that 02 is 
even under m x and odd under m Zl which indicates that 
P has to be even under m x and odd under m z . This can 
only happen if P lies along the c direction, as observed^ 
Finally, for MWO, we see that 0\<j\ is odd under m y . 
This indicates that P 7 also has to be odd under m y . In 
other words P can only be oriented along the b direc- 
tion, again as observed^*! In this connection one should 
note that this conclusion is a result of crystal symme- 
try, assuming that the magnetic structure results from 
two continuous transitions, so that representation the- 
ory is relevant. This conclusion is at variance with the 
argument given by Heyer et al£& who "expect a polar- 
ization in the plane spanned by the easy axis and the 
b axis which they justify on the basis of the spi- 
ral modeliiSii^ It should be noted that their observation 
that the spontaneous polarization has a nonzero com- 
ponent along the a-axis at zero applied magnetic field 
contradicts the symmetry analysis given here. The au- 
thors mention that some of the unexpected behavior they 
observe might possibly be attributed to a small content 
of impurities. 

It is important to realize that the above results are a 
consequence of crystal symmetry. In view of that, it is 
not sensible to claim that the fact that a theory gives the 
result that the polarization lies along b makes it more 
plausible than some competing theory. The point is that 
any model, if analyzed correctly, must give the correct 
orientation for P. 

It is also worth noting that this phenomenology has 
some semiquantitative predictions. To see this, we mini- 
mize P n0 n + Pnt with respect to P to get 

P7 = -XE,~t r -<v>°~< sin(0> - </><) . (147) 

This result indicates that near the magneto-ferroelectric 
phase transition of NVO one has P oc o^aijSl or since 
the high-temperature order parameter 174 is more or less 
saturated when the ferroelectric phase is entered, one has 
P oc <7i, where <j\ is the order parameter of the lower 
temperature incommensurate phase. 

As we discussed, in the low temperature incommensu- 
rate phase one will have arbitrary signs of the two or- 
der parameters. However, the presence of a smll electric 
field will favor one particular sign of the polarization and 
hence, by Eq. (|147fl one particular sign for the product 
a > a < . Presumably this could be tested by a neutron 
diffraction experiment. 

C. TM025 

The case of TM025 is somewhat different. Here we 
have only a single irrep. One expects that as the temper- 
ature is lowered, ordering into an incommensurate state 
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will take place, but the quadratic terms in the free en- 
ergy do not select a direction in <Xi-cr 2 space. At present 
the data has not been analyzed to say which direction 
is favored at temperature just below the highest order- 
ing temperature. (For YM025, as mentioned above, the 
direction ct\ = is favored.) As the temperature is re- 
duced, it is not possible for another representation to 
appear because only one irrep is involved. However, or- 
dering according to a second eigenvalue could occur. Wc 
first analyze the situation assuming that we have only a 
single doubly degenerate eigenvalue. In this case we can 
have a spin distribution [as given in Eq. I|119|l ] involving 
the two order parameters tj\ and <x 2 which measure the 
amplitude and phase of the ordering of the eigenvector of 
the second and third columns of Table IXVll respectively. 
In terms of these order parameters, the magneto-electric 
coupling can be written as 

-Pint = a mni o* n a m P~i , (148) 

where 7 = x,y,z and n,m = 1,2 label the columns of 
the irrep labeled <Ti and <x 2 , respectively, in Table IXVll 
Since reality requires that a nmi = flj,„ 7 , this interaction 
is of the form 



D. CFO 

Again we start with the familiar magneto-electric in- 
teraction 

Pint — ^ A nmi (T n cr* n P 7 , (151) 

where reality implies that A nmi = A* mn ^. Since we have 
Tcr n = cr*, we eliminate terms with n — m: we need 
two irreps for ferroelectricity. Indeed, the higher tem- 
perature phase with a single order parameter &\ is not 
ferroelectric £ Thus the magnetoelectric interaction must 
be of the form 

Pint = ^2[a 1 (T 1 cr* 2 +a*a* 1 (7 2 ]P 1 . (152) 

7 

Inversion symmetry indicates that a 7 = — a* , so we write 

Fi n t = i ^2 r 1 [cTi(T* 2 - <Tt<7 2 ] Pj 

7 

= 2 2_. r 1 a\G2 sin(0 2 — 0i)P 7 , (153) 

7 



ai 7 |cn| 2 + a 27 |<x 2 | 2 



+6 7 <Ti<T 2 + b*(T* 1 <T2 



(149) 



Now use invariance under inversion, taking note of Eq. 
[)120[l . One sees that under inversion cricr^P^ changes 
sign, so the only terms which survive lead to the result 



F h 



5> 7 P 7 [ 

7 



0-1 



k 2 | 2 ] 



(150) 



As we discussed in connection with Eq. (|137fl we have two 
scenarios depending on whether B — \C\ in Eq. (jl36[) is 
positive or negative. If it is positive, then only one order 
parameter is nonzero and we have a nonzero spontaneous 
polarization according to Eq. (|150fl . In that case, using 
Eq. (|120|l we see that [|<Ti| 2 — |er 2 | 2 ] is even under m x and 
odd under m y . For Pi nt to be invariant under inversion 
therefore requires that P 7 be odd under m y and even un- 
der m x , so P has to be along b as is foundS In the other 
scenario, when B — \C\ is negative, then the right-hand 
side of Eq. H150(l is zero and the state is not ferroelectric. 
For TM025 we are probably in the ferroelectric scenario. 
Our analysis therefore suggests that the spin structure of 
TM025 should be given by Eq. [jll9|) with only one of 
the order parameters nonzero. It would be interesting to 
analyze the diffraction data to test this assertion. 



where r 7 is real. Now use Eq. (|129fl which gives that 
ct\ct\ changes sign under 2 X . So for the interaction to 
be invariant under 2 X (as it must be), P 7 has to be odd 
under 2 X . This means that P has to be perpendicular 
to the x axis. Note that symmetry does not force P 
to lie along the three-fold axis because the orientation 
of the incommensurate wavevector has broken the three- 
fold symmetry. 

In the above analysis we did not mention the fact that 
the existence of the ferroelectric phase rquires a magnetic 
field of about 8-10T oriented along the three-fold axis. In 
principle one should expand the free energy in powers of 
H. Then presumably as a function of H one reaches 
a regime where first one incommensurate phase orders 
and then at a lower temperature the second incommensu- 
rate order parameter appears. Then the phenomenology 
of the trilinear magnetoelectric interaction would come 
into play as analyzed above. There is one additional 
point which merits attention. Namely, q could assume 
a symmetry-related value obtained by one or two three- 
fold rotations about the c axis. In the absence of any 
external perturbation to break the three-fold symmetry, 
the system would spontaneously break symmetry by arbi- 
trarily selecting one of the wavevectors in the star. How- 
ever, it is interesting to speculate whether the application 
of a weak in-plane magnetic (or electric) field would be 
enough to enforce the selection of one of the wavevectors 
in the star. If this were so, then the transverse compo- 
nent of the polarization (which, however, might be small) 
would be rotated by the application of such a small ex- 
ternal field. 
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E. High Magnetic Field 

We can also say a word or two about what happens 
when a magnetic field is applied. In TMO, for instance, 
one finds^ that for applied magnetic fields above about 
10T in either the a or b directions, the lower temperature 
incommensurate phase has a spontaneous polarization 
along the a axis. Keep in mind that we want to identify 
this phase with two irreps and from the phase diagram 
we know that the higher temperature incommensurate 
phase is maintained into this high field regime. So the 
higher temperature phase is still that of T 3 at these high 
fields. Referring to Table IXTfl we see that to get CT m cr* to 
be odd under m x and even under m z (in order to get a 
polarization along the a axis) we can only combine irrep 
Ti with the assumed preexisting T3. Therefore it is clear 
that the magnetic structure has to change at the same 
time that direction of spontaneous polarization changes 
as a function of applied magnetic fieldii^ It is also inter- 
esting, in this connection to speculate on what happens 
if the lower additional irrep had been 1?4 so that T4 and 
T3 would coexist. In that case a^a\ is odd under both m x 
and m z . These conditions are not consistent with any di- 
rection of polarization, so in this hypothetical case, even 
though we have two irreps and break inversion symmetry, 
a polar vector (such as the spontaneous polarization) is 
not allowed^ 

For MWO a magnetic field along the b axis of about 
10T causes the spontaneous polarization to switch its di- 
rection from along the b-axis to along the a axis4& We 
have no phenomenological explanation of this behavior at 
present. This behavior seems to imply that the wavevec- 
tor for H > 10T is no longer of the form q = (q x ,^,q z ). 

F. Discussion 

What is to be learned from the symmetry analysis of 
the magnetoelectric interactions? Perhaps the most im- 
portant point to keep in mind is to recognize which re- 
sults are purely a result of crystal symmetry and which 
are model dependent. For instance, as we have seen, the 
direction of the spontaneous polarization is usually a re- 
sult of crystal symmetry. So the fact that a microscopic 
theory leads to the observed direction of the polarization 
does not lend credence to one model as opposed to an- 
other. In a semiquantitative vein, one can say that sym- 
metry alone predicts that near the combined magneto- 
electric phase transition P will be approximately pro- 
portional to the order parameter raised to the nth power, 
where the value of n is a result of symmetry, (n = 1 for 
NVO or TMO, whereas n = 2 for TM025). 

It also goes without saying that our phenomenological 
results are supposed to apply generally, independently 
of what microscopic mechanism might be operative for 
the system in question. (A number of such microscopic 
calculations have appeared recentlyii2i£2iffii£iiS) There- 
fore, we treat YM025 and NVO with the same method- 



ology although these systems are saidiS to have differ- 
ent microscopic mechanisms. A popular phenomenolog- 
ical description is that given by Mostovoyi^ based on 
a continuum formulation. However, this development, 
although appealing in its simplicity, does not correctly 
capture the symmetry of several systems because it com- 
pletely ignores the effect of the different possible sym- 
metries within the magnetic unit cell. 63 Furthermore, it 
does not apply to multiferroic systems, such as TM025 
or YM025, in which the plane of rotation of the spins 
is perpendicular to the wavevectorii (The spin-current 
modeii£ also does not explain ferroelectricity in these 
systems.) In addition, a big advantage of the symme- 
try analysis presented here concerns small perturbations. 
While the structure of NVO and TMO is predominantly 
a spiral in the ferroelectric phase, one can speculate on 
whether there are small spiral-like components in the 
nonfcrroelectric phase. In other words, could small trans- 
verse components lead to a small (maybe too small for 
current experiments to see) spontaneous polarization? If 
we take into account the small magnetic moments in- 
duced on the oxygen ions, could these lead to a small 
spontaneous polarization in an otherwise nonferroelectric 
phase? The answer to these questions is obvious within 
a symmetry analysis like that we have given: these in- 
duced effects still are governed by the symmetry of the 
phase which can only be lowered by a spontaneous sym- 
metry breaking (which we only expect if we cross a phase 
boundary). Therefore all such possible induced effects 
are taken into account by our symmetry analysis. 

Finally, we note that the form of the magneto-electric 
interaction ~ M 2 P suggests a microscopic mechanism 
that has general validity, although it is not necessarily the 
dominant mechanism. This observation stimulated an 
investigation of the spin phonon interaction one obtains 
by considering the exchange Hamilton 

U = J2 J ^(hj)S a (i)S p (j) (154) 

when J a p(i,j) is expanded to linear order in phonon 
displacements^ After some algebra it was shown that 
the results for the direction of the induced sponta- 
neous polarization (when the spins are ordered appro- 
priately) agrees with the results of the symmetry argu- 
ments used here. In addition a first-principles calculation 
of the phonon modes led to plausible guesses as to which 
phonon modes play the key role in the magneto-electric 
coupling. But whatever the microscopic model, the phe- 
nomenology presented here should apply. 

V. CONCLUSION 

In this paper we have shown in detail how one can 
describe the symmetry of magnetic and magneto-electric 
phenomena and have illustrated the technique by dis- 
cussing several examples recently considered in the liter- 
ature. 
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The principal results of this work are 

• We discussed a method alternative to the tradi- 
tional one (called representation analysis) for construct- 
ing allowed spin functions which describe incommensu- 
rate magnetic ordering. In many cases this technique can 
be especially simple and does not require an understand- 
ing of group theory. 

• For systems with a center of inversion symme- 
try, whether the simple method mentioned above or 
the more traditional traditional representation formal- 
ism is used, it is essential to further include the restric- 
tions imposed by inversion symmetry, as we pointed out 
previouslyii^Sii 

• We have illustrated this technique by applying it 
to systematize the magnetic structure analysis of several 
multiferroics many of which had not been analyzed using 
inversion symmetry. 

• By considering several examples of multiferroics we 
further illustrated the general applicability of the trilin- 
ear magneto-electric coupling of the form M(q)M(— q)P, 
where M(q) is the magnetization at wavevector q and P 
is the uniform spontaneous polarization. Usually, when 
expressed in terms of order parameters, this interaction 
requires two active irreducible representations. 

• For TbMn2 05 we analyzed the fourth order terms 
in the Landau expansion and predict that the fact that 
the system magnetically orders into a ferroelectric phase 
indicates that the spin state is described by a single order 
parameter according to Eq. i|119[l . The analysis of the 
diffraction data to test this assertion has not yet been 
done. 

• We briefly discussed the implications of symmetry in 
assessing the role of various models proposed for multi- 
ferroics. 
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APPENDIX A: FORM OF EIGENVECTOR 



In this appendix we show that the matrix G of the 
form of Eq. i|83[) [and this includes as a subcase the form 
of Eq. (|80|) ] has eigenvectors of the form given in Eq. 
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We find that 
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where a' and a" are the real and imaginary parts, re- 
spectively of a and similarly for the other complex vari- 
ables. Note that we have transformed the original matrix 
into a real symmetric matrix. Any eigenvector (which we 
denote \R}) of the transformed matrix has real- valued 
components and thus satisfies the equation 



U -1 GU|ii) = \r\R), 



from which it follows that 



[G}V\R) = \rU\R), 



(A3) 



(A4) 



so that any eigenvector of G is of the form U|i?), where 
all components of \R) are real. If \R) has components 
rl, r2, . . . r7, then 

U\R) = [rl,r2,r3,{r 4 + ir 5 )/y/2,(r4-ir 5 )/y/2, 

{r 6 + ir 7 )/V2, (r 6 - ir 7 )/y/2] , (A5) 

which has the form asserted. 



APPENDIX B: IRREPS FOR TM025 

In this appendix we give the representation analysis 
for TbMn2 05 for wavevectors of the form (i, 0, q), where 
q has a nonspecial value. The operators we consider are 
E, m Xl m y and m x m y , as defined in Table IxTnl Note 
that niy(x, y, z) = (x + 1, y, z), so that = —1 for this 
wavevector. Thus, the above set of four operators do 
not actually form a group. Accordingly we consider the 
double group which follows by introducing — E defined 
by m 2 y = -E, (-E) 2 = E, and (-E)O(-E) = O. Since 
addition has no meaning within a group we do not discuss 
additive properties such as (E) + (-E) = 0. Then, if we 
define —O = (-E)O, we have the character table given 
in Table IXlXl 
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TABLE XIX: Character table for the double group of the 
wavevector. In the first line we list the five classes of operators 
for this group. In the last line we indicate the characters for 
the group G which is induced by the n-dimensional reducible 
representation in the space of the a spin component of spins 
in a given Wyckoff orbit. 



The Mn 4+ Wyckoff orbits contain two atoms and all 
the other orbits contain four atoms. In either case we 
may consider separately an orbit and a single component, 
x, y, or z of spin. So the corresponding spin functions 
form a basis set of n vectors, where n = 2 for the sin- 
gle spin components of Mn 4+ and n = 4, otherwise. In 
each case, the operations involving m x and/or m y inter- 
change sites and therefore have zero diagonal elements. 
Their character, which is their trace within this space of 
n vectors is therefore zero. On the other hand E and 
—E give diagonal elements of +1 and —1, respectively. 
So their character (or trace) is ±n and we have the last 
line of the table for this reducible representation G. 
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In this character table we also list (in the last line) the 
characters of these operations within the vector space of 
wavefunctions of a given spin component over a Wyckoff 
orbit of n sites. Comparing this last line of the table to 
the character of the irreps we see that G contains only the 
irrep T2 and it contains this irrep n/2 times. This means 
that for the system of three spin components over 12 
sites, we have 36 complex components and these function 
generate a reducible representation which contains T2 18 
times. If there were no other symmetries to consider, 
this result would imply that to determine the structure 
one would have to hx the 18 complex-valued parameters. 
The two dimensional representation can be realized by 
Eq. i|121|) ■ The basis vectors which transform as the first 
and second columns, respectively of the two dimensional 
representation are given in Table IX*Xl One can check the 
entries of this table by verifying that the effect of m x and 
m y on the vectors of this table are in conformity with Eq. 



TABLE XX: Spin functions (i. e. unit cell Fourier coeffi- 
cients) determined by standard representation analysis with- 
out invoking symmetry operations that relate q and q. The 
second and third columns give the functions which transform 
according to the first and second column of the two dimen- 
sional irrep. These coefficients are all complex parameters. 



However, after taking account of inversion symmetry 
we have only 18 real-valued structural parameters of Ta- 
ble Q^n] to determine. 
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